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Abstract 

It is generally assumed that a condensate of paired fermions at equilibrium is characterized by a macroscopic wave- 
function with a well-defined, immutable phase. In reality, all systems have a finite size and are prepared at non-zero 
temperature; the condensate has then a finite coherence time, even when the system is isolated in its evolution and 
the particle number N is fixed. The loss of phase memory is due to interactions of the condensate with the excited 
modes that constitute a dephasing environment. This fundamental effect, crucial for applications using the condensate 
of pairs’ macroscopic coherence, was scarcely studied. We link the coherence time to the condensate phase dynamics, 
and we show with a microscopic theory that the time derivative of the condensate phase operator Oq is proportional to 
a chemical potential operator that we construct including both the pair-breaking and pair-motion excitation branches. 
In a single realization of energy E, Oq evolves at long times as h where is the microcanonical chem¬ 

ical potential; energy fluctuations from one realization to the other then lead to a ballistic spreading of the phase and 
to a Gaussian decay of the temporal coherence function with a characteristic time oc In the absence of energy 
fluctuations, the coherence time scales as N due to the diffusive motion of Oq. We propose a method to measure the 
coherence time with ultracold atoms, which we predict to be tens of milliseconds for the canonical ensemble unitary 
Fermi gas. 
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1. Setting the stage 


Coherent gases of ultracold atoms confined in immaterial non-dissipative traps are unique examples of isolated 
macroscopic quantum systems. The value of their intrinsic coherence time is therefore a fundamental question. But 
it is also a practical issue for all the applications which exploit macroscopic coherence, such as interferometry or 
quantum engineering where one generates non-trivial entangled states ^ coherent evolution lilli . Coherence time 
measurements are presently being performed in cold Bose gases |lj,l5ll^. Experiments on Ferm i gas es, which up 
to now focused on traditional aspects of the A-body problem, such as thermodynamic properties are moving 

towards correlation and coherence measurements i9l|. T his turn will open a new research field, including the strong 
coupling regime : that of fermionic quantum optics llOl] . However, a theory predic ting the coherence time of a pair- 
condensed Fermi gas was missing, except in the limiting case of zero temperature 111 ill . In this paper we present the 
first microscopic theory bridging this theoretical gap in a general way. Our approach holds for other physical systems, 
such as mesoscopic Josephson Junctions, provided that the environment-induced decoherence is sufficiently reduced. 

For a Bose-condensed gas of bosons, the finite coherence time is due to the spreading of the probability distribution 
of the condensate phase change during t. At zero temperature and in presence of interactions, a ballistic phase 
spreading is caused by atom number fluctuations in the sample. This effect has been measured by interfering two 
initially mutually-coherent condensates, whose particle number fluctuates due to partition noise Contrarily to 

lasers, which are open quantum systems, and somehow unexpectedly, a ballistic spreading persists in Bose-Einstein 
condensates for a fixed atom number at non-zero temperature lll2lll3l] . Fluctuations of the energy, another conserved 
quantity, then play the same role as number fluctuations. 

For an unpolarized pair-condensed Fermi gas, the study of coherence time presupposes a clear definition of the 
condensate phase and of the corresponding operator 0q ll ill . Furthermore, at non-zero temperature the speed of 
variation of the phase should include the contribution of two excitation branches; the fermionic pair-breaking one and 
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the bosonic one exciting the pair motion. For the fermionic branch Anderson’s Random Phase Approximation (RPA) 
1 141] is enough. For the bosonic branch however, we need the equivalent for fermions of the Bogoliubov method to 
construct quasiparticle creation b\ and annihilation ba operators and to express d0o/df in term of these operators. 
More than that, we need to include interactions among quasiparticles in the evolution of the ba- This is a non-trivial 
consequence of the dependence of condensate wavefunction on the total number of particles N even for a spatially 
homogeneous system, and clearly goes beyond the RPA program. 


2. Correlation function decay 


Below the critical temperature, the time-correlation function of the pairing field i^|(r)i^|^(r') where i^cr(r) is the 
fermionic field operator of the spin cr component, is dominated at long times by the condensate contribution: 

gi(t) = {aliOam) (1) 

where ao - f d^rd^r'^J(r, r')iAi(J’)<AT(*') is ths component of the pairing field on the condensate wavefunction in. At 
equilibrium the system is in a mixture of A-body eigenstates with weights n,i. We therefore study the correlation 
function gj(f) in the eigenstate of energy and particle number N^. To exploit the weak relative fluctuations in 
the number of condensed pairs for a large system, we split oq into modulus and phase hermitian operators 

flo = (2) 

and we approximate the number of condensed pairs operator Ao by its mean value Aq in the equilibrium state to obtain 

gf(f) - (3) 


The operator W, difference between H transformed by e‘®" and H, 


W = e"‘®“Ae‘®» - A = -i[0o, H] - ^[0o, [Bo, H]] + ... 


1 , 


(4) 


is approximatively A times smaller than H. Indeed e‘®“, like ao, changes the total particle number by a quantity 0{N'^). 
While H is an extensive observable, W is intensive and the double commutator in (|4]i is of order 1 /A. In equation Q 
formally appears the evolution operator of the Hamiltonian H perturbed by W, and restricted to the eigenstate of 
H. Up to a phase factor, the function gj /Aq is then proportional to the probability amplitude that the system prepared in 
is still in that state after a time t in the presence of the perturbation W. A standard way to obtain a non-perturbative 
approximation of this amplitude is to use the Green function or the resolvent operator G(z) - (zl - {H + W)) of the 
perturbed Hamiltonian. Within the projectors method (see §III.B.2 of S), we introduce an effective non hermitian 

Hamiltonian Aeff(z) governing the evolution restricted to id est {>!/x\G{z)\4f a) - . This leads 

to 


g\{t) - Ao £ 


dz 




with {H,ff(z))A ^Ea + <W), + {WQa- 


Qa 


-QaW)a (5) 


jc2i7Tz-{H,ff(z))A ' .. ~"zQa-Qa(H+W)Qa' 

The integration domain C in the complex plane is any straight line z - x + ir], rj > 0 fixed, where the real number 
X runs from -i-oo to -oo. The notation (A)^ s {iPa\M>Pa) was used, and the operator Qa - 1 - that projects 

orthogonally to Itf/A), was introduced. Keeping in (Aeff(z)),) terms up to order two in W, hence omitting QaWQ^ in the 
denominator, and neglecting the z dependence (Aefr(z))/i ~ {H^sipA + iO^))i (pole approximation), we obtainfj 13 


glit) 


: with h{5A - ija) = (WQa 


Qa 


(Ea + mQA - QaHQ. 


■QaW}a 


( 6 ) 


1. The pole approximation implicitly assumes that is coupled to a broad energy continuum S. As a consequence the exponential law in 
© holds only at times longer than the inverse frequency width of the continuum, i.e. longer than the quasi particle correlation time introduced 
below. 

2. In order to have a branch cut in the resolvent and a pole in its analytic continuation, which is required to justify rigorously the exis¬ 
tence of a nonzero imaginary part yx and to give a precise meaning to Eq. one must take the thermodynamic limit (A ^ +oo for a 

fixed density). We face here an unusual subtlety: as shown by a generalisation of the reasoning around Eq. {H, the shift function fT3] 
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Remarkably, is the decay rate of the state induced by W as predicted by Fermi’s golden rule. 

The leading term under the exponential in ® is (W)i, it is of order like W. A key step in its interpretation is 
to remark that, according to the expansion in (|4]i, in the Heisenberg picture 

At this stage it may seem difficult to obtain a tractable explicit expression of d0o/df and to go beyond a purely formal 
result for the phase dynamics. Fortunately this is not the case and, as we will show in the next section, the coarse 
grained time average of d0o/df in a weakly excited gas is proportional to a chemical potential operator, which is in 
essence a thermodynamic quantity: 


hd6o 

2'd7 


s=F,B a 


de. 


dN 


( 8 ) 


The sum on the right hand side runs over both the gapped quasiparticles fermionic branch of excitation (in the homo¬ 
geneous case a includes both an orbital and a spin index, a = k, cr) and the bosonic one which, in the thermodynamic 
limit and for an homogeneous system, has a phononic behaviour (a is then only orbital, a - q). By requiring that the 
gas is weakly excited we mean that the thermal depletion of the condensate of pairs must be small. This requires in 
particular that the average number of quasiparticles is a small fraction of the total particle number. The coarse grained 
time average is taken over a time long with respect to the inverse of the quasiparticle eigenfrequencies yet short 
with respect to the typical time-scale of variation of the quasiparticle number operators hs,a, which is allowed if the 
quasiparticles are in the weakly collisional regime. Finally fioiN) is the zero temperature chemical potential of the 
gas with N particles, that is the derivative of the ground state energy with respect to N. We interpret the second term 
on the right hand side of (l8]l as a “chemical potential operator" in the sense that its quantum average is the adiabatic 
derivative of the quasiparticle gas energy Yjs=f,b Ha ^s,a{ns,a) with respect to N, that is at fixed quasiparticle popula¬ 
tions {hs,a)- Equation (O establishes the link between the phase derivative and the chemical potential at the level of 
quantum mechanical operators in a multimode microscopic theory. In that respect, it goes beyond the usual second 
Josephson relation for the phase of the superconducting order parameter (see §3.4 in reference ill]). 

By taking the average of equation ([8]l in the stationary state and using the Eigenstate thermalization hy¬ 
pothesis IB to identify the quantum average in an eigenstate with the microcanonical average, we recognize the 
microcanonical chemical potential jUmc at energy and particle number and obtain: 

= -2;U„e(E,i,A^i) (9) 

at 


A^(z) = (WQ^-x —— ^^~x~Qa^)a^ with Imz > 0, tends to zero as 1/A^, and so does <5^ - in (6l. The way out is to introduce a time 

QAZ+QAiF^-H)Q^ 

scale of order N, setting t = Nt with t fixed at the thermodynamic limit. One performs the change of variable z = Ea + + z'/N in the integral 

of and one chooses t] = t]' fN {?]' > 0 fixed) in the integration domain C = {z = x + ir], x = +oo -oo}, which leads to the integration domain 
C' = {z' = x' + iT]',x' = +00 -oo} over z', and one neglects the terni {W),i in the denominator of the shift function, as well as the terni Rqz'/N 
(to be consistent with the omission of QaWQa). The factor 1/A^ in the element of integration dz = dz'IN allows one to pull out a factor 1/A^ in the 
denominator of the integrand, and one obtains in the thermodynamic limit the following result. 


lim -=- 

A'-^ + CO 


r dz' 

Jc 2k z’ - 


= e with = lim NA^iiri' /N) 

N--*+oo 


which is 77 '-independent due to the analyticity of the integrand in the half-plane Imz' > 0. For a large but finite size system, this leads to the 
definition N{6;i - \y,\) = If one rather takes the thermodynamic limit at fixed t before setting t to infinity, one obtains, after expansion of the 
integrand of s to first order in the shift function, another definition: limAi ^+00 /VlVo- 1 ] ~ i-itmSA - m) with N{Sa - m) = 
lim^^o+ lim/v-^+oo A^A^(i 77 ). To show the equivalence of these two definitions, let us inject in the shift function a closure relation on the eigenstates 
of H, A^(z) = Ea - E^), and let us show that the granularity of the distribution of E^ - Ex in this weighted sum tends to 

zero faster than t]' jN. To this end, we view each eigenstate as a coherent superposition of a central phonon Fock state and of a weak halo of phonon 
Fock states that emanate from the central one by Beliaev-Landau processes 1 phonon ^ 2 phonons. Even if W preserves the phonon number, 
see Eq. it can couple the central Fock state of to the halo of At worst, a single Beliaev-Landau process is involved, in which case 
E^ - Ex = ±(fB,ki + ” ^s,ki-bk->)^ where ki and k 2 are the wavevectors of the two emitted or absorbed phonons; since the pair {ki,k 2 } takes a 

number oc of different values, where V is the volume of the system, E^ - Ex varies by steps oc l/V^, which is indeed «: t}' IN. 
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We omitted here the coarse grained time average as the quantum average is taken in an exact eigenstate of the system. 

The next term under the exponential in ® is of order IjN and subleading. In order to show that, we express this 
term in terms of the correlation function of d0o/df in up to a contribution of order IjN^, 

7dflo(Od0o(O)\ /d0o^"^ 

/ A 


7a + i^A 


f 


dt 




dt dt 


df 


( 10 ) 


This is equivalent to (|6l) as can be checked by inserting a closure relation on the eigenstates of N. The t - 0 value 
of the integrand is Var^(d0o/df) = (9(1/A^) (this comes from adding up the variances of independent quasiparticles 
numbers in the canonical ensemble and it overestimates the microcanonical variance); the function then decays in a 
time Tc which is the typical collision time of quasiparticles and hence the correlation time of the hs^a- Altogether we 
estimate l 7 i + i(5.i| ~ TeVar(d0o/df) = 0(1 IN). The energy shift 6 a is thus of the same ordering as the subleading term 
[^ 0 , [Oq, ^]] in W, i.e. N times smaller than (|9l); we neglect both terms for a large system. In contrast, we keep ja, since 
it is the only term which leads to an exponential decay of the correlation function g'^. Eq. (fTOl i provides a physical 
interpretation of ja, if one remembers that, in the brownian motion theory, the integral of the velocity correlation 
function gives the position diffusion coefficient, ja is thus the phase diffusion coefficient of the condensate of pairs 
when the system is prepared in the microcanonical ensemble corresponding to Ii/za), hence the notation ja - D(Ea,Na). 
We finally keep 


g\{t) 


me 

( 11 ) 

Eq. ([Tol l even gives a way of calculating D: If one can write kinetic equations for the quasiparticles numbers appearing 


in d0o/df, see Eq. ij, one can calculate their time correlation functions as done for bosons in reference IlSll and 
conclude that 

D(E,N)^-(pX)- (12) 
The matrix M, with coefficients Msa.s’a’, is the matrix of the linearised kinetic equations ^6n - M6n that give the 
evolution of the fluctuations 6ns,a of the quasiparticle numbers, collected in a single vector 6ii, around their stationary 
values hs^a- To define the other notations, one introduces as in MIST the covariance matrix Ccan of the quasiparticle 
numbers n , q. in the canonical ensemble with N particles and a mean energy E, the energy vector e with components es^a 
and its dual vector e oc Ccan? normalised such that e-e - 1.0 The vector A - (2 IK) del dN then collects the coefficients 
of hs^a in the expression (l8]l of -d0o/df, the matrix P is the non-orthogonal projector^ such that P6n - 6n — e(e ■ 6n) 
for all 6n, and the microcanical covariance matrix is Cmc = P^C^anP Ei. 

The final step is to take the statistical average of Eq. (fTTT i over the probability distribution 11^ of the states that 
constitute the mixed state of the system. For large N, we assume that energy and atom number fluctuations around the 
mean values E and N are weak in relative value. This is the case if, for example, describes a canonical or grand 
canonical ensemble. We assume Gaussian fluctuations and linearize /ime around (E,N) while, to this order, we keep 
only the central value D(E, N) of the next-to-leading term. We are led to the calculation of a Gaussian integral with a 
phase factor exp{2i[(diJ.mciE, N)IdE)(E - E) + (dp.mc(E, N)!dN)(N - N)]t!h). Altogether this leads to the main result 
of this work: 


gi(t) No& 






(13) 


In presence of energy or atom number fluctuations, the thermal blurring at long times consists in a Gaussian decay of 
the correlation function giit), with a characteristic time 


(2fbr//i)"" = Var A 




™(£,A)h-£^^' 


me / -Tj 


dN 


dE 


{E,N) 


(14) 


which diverges as for normal fluctuations. The phase diffusion coefficient D leads to an exponential decay with a 
characteristic time diverging as A. As expected it is a subleading effect at long times, except if the system is prepared 
in the microcanonical ensemble in which case the intrinsic phase diffusion may be directly observed. 


3. One has (Ccan).vQ', 5 '(i' = ± ^s,a) and hs,a = 1 T i] where the upper (lower) sign holds for the bosonic (fermionic) 

excitation branch. 6 and ^ are left and right eigenvectors of M with a zero eigenvalue, as shown by the conservation of energy and the inspiring 
rewriting ^ oc dn/dT, where T is the canonical ensemble temperature and n is the vector with components hs^a B. 

4. Pe = (1 and projects onto the subspace of fluctuations 6n of zero energy, \6n | e ■ dn = 0|, inside which one defines the inverse of the matrix 
M. 
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3. Microscopic derivation of the phase operator equation 


We give here the first (to our knowledge) microscopic derivation of equation ([8]l, relating the evolution of the 
phase operator of a pair-condensed gas to a “chemical potential operator". 

The contribution of the fermionic branch of excitations to d0o/df can be obtained from linearized equations of 
motion for small fluctuations of the pair operators lAjlAj and around the mean-field state in Anderson’s 

Random Phase Approximation (RPA) 11411 . Using equation (120) of reference Q to extract the time average of 
d0o/df, and rewriting equation (86) of the same reference in terms of the fermionic quasiparticle occupation numbers 
hf^a, we get! 


h d0o 
2'd7 


= ^(N) + %(N- 

RPA dN 




t>'=k,a' 






(15) 


where fF.k.o- is the BCS excitation spectrum of an homogeneous system and N the BCS average particle number in the 
grand canonical ensemble of chemical potential /r.Q 

We encountered fundamental difficulties in deriving the phonon branch contribution to equation (O within the 
RPA.0 We therefore decided to treat the problem variationally with the most general time-dependent pair coherent 
state Ansatz:^ 


|i/r) = Af(f)exp 


^r(r, r';f)iAj(r)<A|(r') 


10 ), 


(16) 


r.r' 


/ 


Here N ensures normalization and the r(r, r') form a set of independent variables. The space has been discretized on 
a cubic lattice of step I, which we take to zero in the end of the calculations. The held operators obey anticommutation 
relations of the kind: {i^o-(r), {j/'^{r')} - 6a-,a-'Sr.r'■ Section §9.9b of reference ifl^ constructs from r(r, r'), r*(r, r') 
the set of canonically conjugate variables <l)(r, r'), ®*(r, r').0 This held O should not be confused with the usual 
pairing held (lAiiAT) LJ When the pairs in (fThl l are at rest, r(r, r', t) depends only on r - r' and the Fourier transform 
of O reproduce the 14 amplitude of the k t, -k J. pair of BCS theory 111 ill , while the Fourier transform of the pairing 


5. We use ^~~ik = ^F.k.t + where Aq and 4 refer to notations of (mi , and we use Eq. (74) of that reference to recognize de/rk o-Zd/r. 

6. To obtain Eq. fl5) in this form we reintroduced the trivial phase term -2fi(N)t/h that does not appear in the reference (m due to the use of 
the grand canonical Hamiltonian. 

7. The RPA result il5\ . restricted to the linear order in the pair operators, does not include the contribution of the phonon branch. One might 

hope to obtain this contribution by pushing the RPA calculation to the quadratic order in the pair operators as follows: First, one computes dOo/dt 
up to the quadratic order. Second, one collects all the RPA pair operators inducing a center of mass momentum change hq, that is «-k-qj.^kT’ 
^k+qT^-ki’ ^k-fqT^*^^^ where ntT annihilates a fermion of wave vector k and spin t, and one writes the matrices Xq of their 

closed linear RPA equations of motion. Each RPA matrix Xq has two collective modes, with opposite energies linear in q for small q. The RPA 
operators are then expanded over the eigenmodes of Xq, with amplitudes b-q and on the collective modes that annihilate and create a collective 
excitation of momentum Thq. Third, one inserts this modal expansion in the RPA operators appearing in dOo/dt and obtains a quadratic expression 
in the modal amplitudes, hence terms in Pqbq as in m- Unfortunately, this RPA approach is not reliable and must be abandoned because the 
RPA operators, although linearly independent, are not quadratically independent, as one can see by rearranging the pair operators using fermionic 
anticommutation relations, e.g. 

("LqT^kTXalk-qi^-kf) = («LqT“-k-qi)(«-ki«kT) 

This shows that (/) there is no unique way of expressing dOo/dt as a quadratic function of the RPA operators, (ii) the modal amplitudes are not 
quadratically independent, (in) the coefficient of is not uniquely determined by this RPA approach. 

8. We use here for simplicity an Ansatz of the time-dependent BCS type in which the number of particles has quantum fluctuations (not to be 

confused with the thermal fluctuations of the grand canonical ensemble). The use of an Ansatz with a flxed number of particles, possible although 
more difficult fT3] . would change the coefficients in the energy functional by a relative coiTection and would not change the spectrum in 

the thermodynamic limit. 

9. Our variables T and ct correspond to the variables z and /3 of reference fl^ when one replaces the quasiparticle operators and vacuum in 
Eq. (9.132) of (H by the particle ones. If E and ^ ai‘e the matrices of elements Pr(r, r') and l^O(r, r') respectively, then ^ = -E(l + E^E) 

according to Eq. (9.144) of fl^ . Their variational equations of motion follow from the usual Lagrangian L = ih [{il/\(d/dt)\i//) - c.c.] /2-(^|7/|^). d> 
is cleverly defined such that L = ihf' Xr,r' r0i9rfl>(r, U) - c.c.] /2-'7i(d>, ffi*), leading to dt^ = (\hl^)~^. This shows that the conjugate 

variable of <I) in the Hamiltonian formalism is for the usual Poisson brackets, that is O* for the Poisson brackets {<!), <!>*} = (i^/^)“^ knowing 

that one must have dt<b = 

10. With the same matrix notation as in note 9 one has = ^(1 - according to Eq. (9.146) of the reference fT^ . 
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field is -t/kW = -^(1 - For moving pairs we have no physical interpretation, but the squared norm of O 

is still half the mean number of particles in the state |i/r): 

I = ||(b||2 = /®2|(l)(r,r';0^ (17) 

r,r' 

The classical Hamiltonian governing the evolution of the field O is given by 

= ( 18 ) 

and may be expressed explicitly using Wick’s theorem. In the following however, we will only need the invariance of 
"K under a global phase change <t)(r, r') ^ e'^<t>(r, r'), Vy6lR(f/(l) symmetry), consequence of the conservation of 
the particle number N by evolution with H. At zero temperature and for a fixed N the field <t)(r, r') is frozen, up to a 
global phase factor, into the minimizer <I>o(r, r') = {NI2y^^(/)Q(r, r') of 'Tf. 0o is chosen real and normalized to one. It 
depends on N even in the spatially homogeneous case and differs from the condensate wavefunction ipo in the same 
way that <t) differs from the pairing field (see note 10). At sufficiently low temperature one can expand 'TT in 

powers of the small deviations of <1) away from the circle y i-> e'^<l>o(r, r'), locus of the minima of "TY for fixed N. To 
this end, we split the field into its components parallel and orthogonal to 0o : 

0(r, r') = e'®[nF20o(r, r') + A(r, r')] (19) 


The phase 6 can reach arbitrarily large values while A is bounded. This framework allows us to develop a sys¬ 
tematic perturbation theory in powers of the field A (see appendix A), the fermionic equivalent of the Bogoliubov 
t/(l)-symmetry conserving approach for bosons 1^ . Provided that A stays small, the phase 6 remains close to the 
condensate phase 9o as we shall see. We therefore write down the equations of motion of 6 and of the fields A, A*. 
At the end of the calculations we systematically eliminate the condensate variables with the relation n - Njl— ||A|p, 
consequence of (fTTl i. and we restrict ourselves to order 2 in A, A*. 

The main challenge of the calculation is the occurrence of a term linear in A, A* in dfl/df, resulting from the 
fact that 00 depends on the number of pairs d. Without this term, one would simply expand the field A on the 
eigenmodes of its small linear oscillations obtained from a quadratization of the Hamiltonian 'H at fixed N: 


/A(r,r';f)\ 

iA*(r,r';f)/ 




lua{r, r')\ 
l^v„(r, r')/ 




A'a(r,r')'i 

l<(r,r')j 


( 20 ) 


where the sum runs over the eigenmodes of positive energy ea, normalized as \\ua\?' — I|vq.|P = 1. To this order, 
bait) - b„(0)e“"^''7fi One would insert the expansion (l20l i into dOjdt and take a coarse grained temporal average to 
get rid of the oscillating terms, dfl/df would then contain the expected linear combination of the numbers of bosonic 
quasiparticles ns^a - In reality, the problem is more subtle: due to the interaction among the quasiparticles, ba 
does not vanish and is of order two in A and A*. The contribution of the linear term in d6>/df is then comparable to 
that of the quadratic ones. It is calculated in the appendix A, using in particular the bounded nature of the field A 
(a consequence of the t/(l)-symmetry preserving nature of the expansion ( fT9l l) and the Hellmann-Feynman theorem. 
One finds 

~ 2 ^ = Fo(A) + 2] + 0(||A||3) (21) 

a 


We now briefly discuss the form of the energy spectrum for a spatially homogeneous system, in the continuous 
limit I 0 for a s-wave contact interaction with a fixed scattering length between opposite spin fermions. For each 
value of the total wave vector q, there exists (/) at most one discrete value eB,q, (U) a continuum parametrized by two 
wave vectors (ki, t; k2, i) i--» e^^ki.T + fF.k,,! of constant sum (ki -t- k2 = q), where e^k.o- is the BCS dispersion relation. 
The branch e^ q coincides with that of reference 112111 , as we have checked. It has a phononic start and corresponds 
to the bosonic elementary excitations of the Fermi gas, whose contribution to the phase dynamics was missing. The 
continuum corresponds to the excitation of two fermionic quasiparticles. Indeed, since the Hamiltonian H contains 
an even number of factors 0 and 0^, each annihilating or creating one quasiparticle, fermionic quasiparticles can only 
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Figure 1: Some experimental considerations: (a) Trapping configuration proposed here to measure gi (t) via Ramsey interferometry: the condensed 
paired Fermi gas is confined in the main trap (with a flat bottom on the figure); one transfers on average at most two atoms (in the form of a dimer) 
in the (very narrow) secondary trap via a resonant tunneling effect, which can be tuned by a barrier of adjustable height; in this way, one creates a 
phase reference, which is made to interfere with the condensate after an evolution time t. (b) Thermal blumng time of a coherent Fermi gas in the 
unitary limit in the canonical ensemble, as a function of temperature T in units of the Fermi temperature Tf = Sflkg. Discs: from the equation of 
state measured in reference d. Dashed line: expression i22\ deduced from an approximated equation of state (see text). 


be created by pairs from the ground state. The corresponding biexcitations are not physically independent.^!] and 
duplicate the contribution of the RPA to d0o/df. They must therefore not be included in Eq. (1211 1. 

Two more remarks are needed to obtain ([^i. (/) The fields O differ and so do the phases 9q and 9. Their 

coarse grained temporal averages, however, only differ by a term of order ||A|p, which, bounded hence negligible in 
the long time limit, does not contribute to the phase blurring of the condensate of pairs. [3 (//) The phase 9 of our 
variational approach is a classical variable, whereas 9q in (|8ll is a quantum operator. This gap can be bridged by 
using the quantization procedure exposed in Chapter 11 of reference lIlQIl where the ba of the bosonic branch are in 
the end replaced by bosonic operators]*^ b„. [ba,b]j\ - 1. We argue that equation (O, linking d0o/df to the chemical 
potential operator, and the resulting equation (fT3l l should hold beyond the validity range of the microscopic variational 
derivation presented above, and should apply even to the strongly interacting regime, provided that the temperature is 
low enough for the quasiparticles lifetime to be much longer than the inverse of their eigenfrequency. Indeed, in the 
limiting case where one can neglect the fermionic excitation branch and drop the non-phononic part of the bosonic 
branch, Eq. (fSll can be derived from the irrotational version of the quantum hydrodynamic theory of Landau and 
Khalatnikov 12211 (see appendix B). 


4. Towards an experimental observation 


Let us briefly explain how an experimental evidence of the thermal blurring of a condensate of pairs could be 
obtained. The key idea is to bosonize the atomic Cooper pairs into deeply bound weakly interacting dimers during the 
preparation and the measurement stage. This can be done in an adiabatic reversible way 02311 by tuning the scattering 
length to a small and positive value thanks to a magnetic Eeshbach resonance. It allows one to (/) produce a sample of 
dimers with weak number fluctuations from a melted Mott phase of an experimental realization of the Bose Hubbard 
model ll^ . {ii) control tunneling between the main trap (containing the N particles) and a very narrow secondary 
trap by adjusting the height of a potential barrier il (FigOi), (///) detect by fluorescence a single dimer in this 
secondary trap. For the measurement of the gi(0 function, we adapt to the case of paired fermions the interferometric 


11. Exciting a = (ki, f;k 2 ,i) and a' = (kj,t;k^, i) amounts to exciting a" = (kj,f; k 2 , i) and a"' = (ki,t;k^, i). 

12. Expressing in {qq) the pairing field in terms of <I), one realizes that, for small A, {oq) = V^(l + D(||A||)) and, since A = 0(||A||^), one 
has Tq + D(||A|p). 

13. More precisely, these operators are bosonic only for a weak density of excitations. For a spatially homogeneous system and in a U{\) 

symmetry breaking formalism (6 = 0), we obtain from Eq.(11.81c) of fl^ extended to the paired case replaced by the BCS fermionic 

quasipai'ticle annihilation operators bko-) and restricted to weakly excited bosonic images (BB^ negligible): = Zk ^k^-ki^k+qT + ^k^k-qT^-ki‘ 

The real coefficients and are linear combinations of the con'esponding Wq and Vq in Eq. and inherit the normalization condition Zk(-^k^^“ 
(y^)2 = 1. Then in a state <16t with 0 = 0, Sq = [&q,^q] - 1 has a mean value and a vai'iance 0(||A|p/A) since ^ko-|iA) - if 

||A||/iV‘/2^0. 
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Ramsey method of reference lE^I . where two Rabi pulses are applied at a time interval t. The bosonized pairs are 
prepared initially in the main trap. A first pulse of angle e transfers on average less than one dimer to the secondary 
trap; in this way, the thermal blurring is not masked by partition noise. Then the system evolves during a time t with 
interactions set to the value at which phase dynamics is to be studied. Last, the gas is rebosonized and a second pulse 
of angle e closes the interferometer before the number «sec of dimers in the secondary trap is measured. The average 
of Usee over the realizations is an oscillating function of the time f, of angular frequency 2/h times the difference of 
the two trapping zones chemical potentials, and of contrast equal to |.gi(0/.?i(0)|. 


5. Explicit results for the unitary gas 


We first estimate the blurring time for a unitary Fermi gas prepared in the canonical ensemble, that is with energy 
fluctuations of variance VarE = ksT^drE. From the equation of state of the unpolarized unitary gas measured in 
reference |@1, and for a spatially homogeneous system (in a flat bottom potential |27]) we find the thermal blurring 
times fbr plotted as discs in figure [T}?. For example, at a temperature T - 0.12 Tf ^ 0.7 Tc, where Tc is the pair 
condensation temperature, we find fbr ~ corresponding to 20 milliseconds for a typical Fermi temperature 

Tf - Bflks - ljuK and a typical atom number A = 10^. As in reference ll2^ . one can also estimate the equation of 
state of the unitary gas from simple dispersion relations for the elementary excitations. For the bosonic branch one 
takes lE^ es.q - with c the T - 0 sound velocity, mc^ — and ^ the Bertsch parameter. For the fermionic 
branch, one takes lE^ = A + - eo)^/(2/o), where A is the gap, and eq and /o give the location of the 

minimum and the curvature of the dispersion relation. Keeping each branch contribution to the mean volumic energy 
to its leading order at low temperature Ef/V ^ and Eg/V ^ ^(kRT)'^ l(hc)^ as in E3l and 

using the experimental values ifiH ^ = 0.376, A = 0.44ef, eq - 0 . 85 ef and the theoretical value Ell /o = 0.846e/r, 
we find 


Nh^ 

(tbrSp)'^ 
\9/2 


0 


(1 +2rr 

(1+0 


( 22 ) 


with 0 = TITf and r= drEfldrEB - e j-giafive weight of the two excitation branches. This 

formula.r^ plotted as a dashed line in figure[TJ), is an exact equivalent to fbr for 0^0. 

At finite times[3 t - O(N^) fbr, the contribution of D to gi(f) in equation (fT3l l is a priori comparable to that of 
since both scale as 1/A in the canonical ensemble. D can be calculated at very low temperatures keeping only the 
phononic part of the spectrum. In the state of the art, it is predicted by various approaches that the bosonic branch is 
convex at low q for the unitary gas 1 3J, lElIll, 


yh?q^ 

fB.q ^ fu:q+ + 0{q^) with y > 0 

9->o am-^c 


(23) 


so that phonons interact through Landau-Beliaev processes 2 phonons <-> 1 phonon as in the weakly interacting Bose 
gas. 0 One can then adapt0 the low-temperature reduction of the general expression (fT^ done in reference ifTll . 


14. The good agreement with the experimental data has to be taken cautiously. If one treats the two branches to all order in kgT, one gets an 

upward shift of more or less constant and equal to 5 over the temperature range of figure 

15. Still the times that we consider are subjected to the constraint t > Tc (see note 1). For the unitary gas Tc-{0.149/©)^ when 0^0, according 
to noteflT] 

16. In the concave case, the leading processes ai'e the scattering events 2 phonons ^ 2 phonons as in reference (23] . 

17. The universalising trick at low temperature is to use, whatever the physical system, mc/h as the unit of wavevector and mc^ as the unit of 
energy. Let us then start from the phase diffusion coefficient and the quasipaiticle congelation time for weakly interacting bosons fl^ . hND!(mc^) ~ 
ciikgT/mc^)^ and mc^rdh ~ C2,p[hl{'\l2mc)\^{mc- jkgT)^, where c\ ^ 0.3036 and C 3 0.05472, and let us review all the coiTective factors 
connecting the weakly interacting Bose gas to the unitary Fermi gas, knowing that D is a quadratic function of PA and a linear function of M~^, see 

Eq. jl2> : (/) for a condensate of pairs, there is an additional factor two in the coefficients A of -dO^ldt , thus a factor four on Z), {ii) the equation of 
state poip) oc p is replaced by p^ip) oc p^^^, where p is the density, so the value of {Nlc)dcldN in A jumps from 1 /2 to 1/3, hence a factor 2/3 on A 
and a factor 4/9 on D, {Hi) PA, being exactly zero for a linear-in-^ bosonic excitation branch due to A oc is proportional at low temperature to the 
dimensionless curvature parameter y, equal to one in [H , hence an extra factor y- in D, (fv) according to quantum hydrodynamics (23] . the reduced 
Beliaev-Landau coupling applitues among the quasiparticles, at fixed values of the reduced wavevectors, are system-dependent only via a global 
factor 1 + A, with A = pPq(p)/[3pq(p)] jumping from 0 to Ay = -1/9, hence a factor (1 + Ay)“^ on the matrix M~^, on D and on the congelation 
time Tc induced by the collisions among the quasiparticles. As a consequence, for the unitary gas, hDNI{mc^) ~ ci[( 47 / 3 )^/(l + Ai,)^]{kBTjme^)^ 
and mc^Tclh~ 03(1 + Al,)^^p[H/{^/2mc)]^{mc^/ksT)^. 
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Furthermore, from a numerical solution of the RPA equations of reference 112 111 we find that y a; 0.1, in agreement 
with |33, 34], so that 


hND 


C0^ 


with 


C ^ 0.4 


ep 0^0 

With this we reach a complete picture of the thermal blurring of the unitary Fermi gas at low temperature. 


(24) 


6. Conclusion 

We have presented the first microscopic theory of the thermal blurring of the phase of a condensate of pairs of 
fermions O, revealing a ballistic blurring and a subleading phase diffusion. The blurring time depends on the 
variance of the total energy of the gas, and on the derivative of the microcanonical chemical potential with respect to 
the energy. The phase diffusion coefficient can be deduced from kinetic equations on the quasiparticle numbers. To 
obtain these results, we have used the fact that the time derivative of the condensate phase is given by the chemical 
potential operator of the gas, see equation We have derived this central relation in a fully microscopic way, 
including both the bosonic and the fermionic branches of excitation. Last, we have proposed a realistic experimental 
protocol to measure this blurring time, that we estimated to be tens of milliseconds for a coherent gas prepared in the 
unitary limit in the canonical ensemble. 
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Appendix A. More on the variational calculation 

Here we derive equation (1211 1 of the time average of 6 within the microscopic model based on the Ansatz (flSI l. 
coherent state of moving pairs. In a first stage one should perform the expansion of the Hamiltonian treating the 
real quantity n and the complex field A in ( fT9l l as independent variables, that is, not fixing the value of ||<t)||. To include 
interactions among the quasiparticles, one must go to third order in A and A*: 

3 

•K(cl), CD*) = To[u, .^o(A)] + ^ Tj[n, 0o(A)](A, A*) + 0(||A|h (25) 

M 

where the tensor Tj is of rank j so that Tj{A, A*) is exactly of order j in A and A*. It may be expressed in terms of the 
differential of order j of '7/ taken at (<1), <!>*) = and restricted to the subspace orthogonal to {<po, 0) and 

(0, (po), with 00 assumed to be real. It does not depend on the phase 0 due to t/(l) symmetry. Eor afixed total number of 
particles, the energy does not vary to first order around the minimizer so that 7“i [A/2,0o(A)] = 0, which is the famous 
gap equation when the system is spatially homogeneous. Eurthermore, one can check that 5„7o[A/2,0o(A)] = 2poiN) 
where poiN) - d£’o(A)/dA is the gas chemical potential at zero temperature, EoiN) - 7o[A/2,0o(A)] being the 
ground state energy. 

The phase and the modulus square of the component of the field ch on the mode 0o are canonically conjugate 
variables, so that -MOfdt = (9„'7T(<1), <!>*). Once this derivative is taken in (l25l l for fixed A and A*, one can fix the 
norm of d) to the value (A/2)'^^ (that is the total particle number is fixed to N), and eliminate n through the identity 
n - ||<l>|p - ||A|p; the field A then remains the only dynamical variable of the problem. The resulting expression is 
useful up to order 2 in A, A*: 

d0 ^ 

“ " 5«^o[A/2, 00(A)] - \\AfdlTo[N/2, 0o(A)] + ^ (l„Tj[A/2,0o(A)](A, A*) + C>(|]A|]3) (26) 

J=i 

The Hamiltonian that determines the evolution of A at fixed particle number is obtained by replacing n with 
N 12— |]A|]^ in (i25l l and by expanding the resulting expression up to order three in A, A*: 

'Kw(A, A*) = £o(A) +f 2 [A](A, A*) +T3[A/2,0o(A)](A, A*) - |]A|]25„T'i[A/2,0o(A)](A, A*) + 0(|]A|]4) (27) 
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with the quadratic form T 2 [A^](A, A*) obtained by subtracting 2^o(A^)l|A|p from Ti[NI2, A*). To compute 

the coarse grained time average of dnT\{Njl, ^q(N)\{A, A*), that is the term linear in A, A* which is problematic in 
d0/df (see the main text of the paper), we write the temporal derivative of the imaginary part of the component of the 
field A on the function (A/2)*^^d0o/dA, 


y4^(f)'''^^(A(r.rVA-(r,n) 


2i 


r.r' 


dA 


Since ifi5,A = I °5a*‘Ha^(A, A*), one gets 
dF 

- 2h— = D ■ -KaCA, a*) = D ■ T2[N](A, A*) + D- -K^^A, A*) + (9(||A|h 
at 

where is the component of TYat of order three in A, A*. We have introduced the differential operator 


® “ Zj i 2 j dA '^A-(r,r')) 

r.r' 


(28) 

(29) 

(30) 

(31) 

(32) 


We shall now take advantage of two identities that exactly hold for all A orthogonal to 
2D ■ f 2 [A](A, A*)=-5„T'i [A/2,0o(A)](A, A*) 

2D ■ 'H'itiA, A*)=2 Af2[At](A, A*) - 5„T2[A/2,0o(A^)](A, A*) + ||A||2527-o[A/2, <^o(A)] 

To prove these relations, one formally considers a field <h' with A + 6N particles and determines in two different ways 
the quadratic expansion of "HN+miA', A'*) in powers of A' and A'*, where A' is up to a global phase as in ([T9l l (written 
for A + dN particles) the component of O' orthogonal to (po{N + 6N). First, one simply replaces A with A + 6N and A 
with A' in dZTl i. and then expands to first order in 6N. The tensor d‘72/dA[A] naturally appears from this expansion. 
Second, one applies to ”7^(0', O'*) the expansion (l25l l around Oo(A) and takes into account the fact that, to first order 
in dA, the component of O' orthogonal to <po{N) contains, in addition to A(l coming from A', a contribution coming 
from dcpoldN: 


A(r,r') = dA^y) ' ^(r,r') + Al(r,r') + (9(dA^ ||A'||^dA||A'||2) 


(33) 


This infinitesimal shift proportional to dA along the direction of dcp^jdA is responsible for the occurrence of the opera¬ 
tor D. Equations OTT i and (l32l i are finally obtained by identification of the two resulting expressions of ’HNi-miA', A'*) 
respectively to first and second order in A' and A'*. 

It remains to combine equations (l26l l. (l29l l. OTT i. and (l3^ after a coarse grained temporal average (over a time 
scale much longer than the oscillation period of the modal amplitudes ba but much shorter than the evolution time of 
the quasiparticle numbers \ba\^), to obtain 


fit =^o(A).^[A](A,A*) 


0(||A|p) 


(34) 


where we used the crucial property that dA/df vanishes (and so does dT/df), since the range of variation of the field 
A is bounded. The quadratic form “7^2 [A] is represented by the matrix rjXXN] with tj - ^ 


and the hermitian scalar product (,) generating the norm || ||: 

f2[A](A,A*) = i/(A,A*),/7X[A] 


0 


-1 


, using a block notation 


(35) 


Then one inserts the modal decomposition (l20l i in the derivative with respect to A of the equation dTSl i: we recall that 
ea and («„, Va), -ea and (v*, «*) are the eigenvalues and the eigenvectors of X.- The coarse grained temporal average 


removes the crossed terms, and the Hellmann-Feynman theorem ensures that0 {(Ua,-Va), 
and finally leads to equation (l2Tl i. 


dN 


) = de„/dA, 


18 . {Ua, -Va) is the dual vector of {Ua, Va) because = V-C- 
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Appendix B. Irrotational quantum hydrodynamics 


In this supplementary section we give an additional derivation of equation (l8]l based on quantum hydrodynamics. 
Although this derivation is not microscopic contrarily to the one presented in the main text of the paper and in appendix 
A, and although it neglects the internal fermionic degrees of freedom treating the pairs of fermions at large spatial 
scales as a bosonic field, it has the advantage of relying only on the equation of state and thus of being valid in all 
interaction regimes. 

To calculate the viscosity of superfluid helium at low temperature, Landau and Khalatnikov have developed in 
1949 the theory of quantum hydrodynamics 12211 . It allows one to determine, to leading order in T, the effect of a 
non-zero temperature on the quantum fluid, at least on the observables that only involve low energy scales and large 
length scales. Remarkably, the only specific ingredient is the zero-temperature equation of state of the fluid, which is 
here the energy per unit volume eoip) of the ground state of the spatially homogeneous system of density p. 

To obtain the time derivative of the phase operator of the condensate of fermion pairs, we refine the theory in two 
ways: by regularizing ultraviolet divergences and by specializing to the irrotational case. 

First, we solve the issue of the Landau-Khalatnikov Hamiltonian ground state energy, that diverges due to the 
zero-point motion of the system eigenmodes. We discretize the space in a cubic lattice of spacing /, a sub-multiple 
of the size L of the quantization volume, which is much smaller than the typical wavelength 2nlq^a of the thermal 
excitations of the fluid but much larger than the mean interparticle distance p 


p 1 <s; 


(36) 


both conditions being compatible at sufficiently low temperature. This is in the spirit of the validity range of hydrody¬ 
namics, which relies on a spatial coarse graining, and it provides a natural cut-off for the wave vectors q by restricting 
them to the first Brillouin zone[3 £) = [-tt//, 7r//[^. In the Hamiltonian one must then replace the differential op¬ 
erators such as the gradient, the divergence and the Laplacian, by their discrete versions, as we shall implicitly do 
below, and introduce the bare energy density eo,o(p), which depends on the lattice spacing 1. Following the ideas of 
renormalization, the zero-point energy of the modes, that formally diverges when I —> 0, adds up to eo,o(p) to exactly 
reconstruct the effective or true energy density eoip), that does not depend on I and is what is measured experimentally. 

Second, we specialize the theory to the case of an irrotational velocity held operator v(r, f) that can then be written 
as the gradient of the phase held operator ^(r, f), itself canonically conjugate to the density field operator p(r, t): 


v(r, t) - —grad ^(r, t) 
m 


with 


[p(r,f),^(r',f)] 


(37) 


This amounts to neglecting the transverse component of the held v(r, t), as done in reference Il22ll to determine the 
phonon-phonon interaction and go forward in the computation of viscosity. In the part icular scale invariant case of 
the unitary Fermi gas, this was justified within the effective held theory in reference 13511 . We note en passant that the 
density of fermionic quasiparticles is exponentially small in l/T and is directly omitted by hydrodynamic theory. 

The steps that follow are rather usual. One starts form the equations of motion of the fields in Heisenberg picture, 
that is the quantum continuity equation and the quantum Euler equation for the potential (whose gradient gives the 
quantum Euler equation for the velocity): 


dtp + div 




= 0 


1 2 

hd,<p^--m\ -po,o(P) 
where {A,B} = AB + BA is the anticommutator of two operators and 

Mo,o(p) = ^eo,o(p) = eoofp) 


(38) 

(39) 


(40) 


19. We also require that the elementary excitations of the system remain phononic over the whole domain D. This imposes hc/I < A so 
mclih > mc^lA » 1 in the BCS limit, and only mcllh > 1 elsewhere, where c is the sound velocity and A the gap at T = 0. The necessary 
condition mcllh > 1 is weaker than the already assumed one, » 1, in the strongly interacting regime or in the BCS limit, since c is then of the 
order of the Fermi velocity. It is more stringent and must be added explicitly in the so-called BEC limit, where the fermion pairs can be considered 
as bosons and form a weakly interacting condensate. 
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is the bare ground state chemical potential at density p. These equations originate from the Hamiltonian 


^hydro “ ^ ^ ] 


• pv + eo,o(p) 


( 41 ) 


The quantum spatial density and phase fluctuations are weak provided that p'^^/ is large enough; the thermal ones are 
weak if in addition q^iJ. is small enough.|3 Under these conditions one can linearize as in 122] the equations of motion 
around the spatially uniform solution: 


p(r,f) = po + 5p(r,f) (42) 

0(r,f) = 0o(f) + 50(rU) (43) 

The operator po reduces to N jL?, where N is the operator giving the total number of particles, and is a constant of 
motion. The operator 0o is the phase operator of the condensate; one has here 

00 = 0o/2 (44) 


since the phase operator 6q in equation (|2]i takes the pairs as the building block, whereas equations (138139b are build 
on the fermionic particles. The spatial fluctuations 6p and b0, of vanishing (discrete) integral over the whole space, 
can be expanded on the plane waves of non-zero wave vector q, and commute with po. One solves the linearized 
equations for bp and b0 and one can use the usual expansion on eigenmodes: 


.,1/2 

bp(r,0-^ Z 


hq 


qe^Z3*n£) 


2otco,o 


1/2 


t 'ipiq r 


(Z?q + B_q)e 


b0(r, f) - 




Po 


t3/2 Au \ 2 tui I ^ 


qefZ^*r)D 


(45) 

(46) 


where the creation operator Bq and the annihilation operator Bq of a phonon with wave vector q and energy hqco,o obey 
bosonic commutation relations [Bq, Bq,] = bq ^- and where we introduced the zero-temperature bare sound velocity 
operator 


Co.o = 


PoPo,o(Po)\'^^ 
m I 


(47) 


It remains to expand the right-hand side of (l39l l up to second order in bp and b0, to extract the zero wave vector Fourier 
component, to perform a coarse grained temporal average to get rid of the oscillating crossed terms BqB_q and BiqBq, 
and to use the identity 


d 

dp 


[PBo,o(P)] 


1/2 _ 


Po,q(p) +PPofi(p) 

2[pp;o(p)]F2 


(48) 


to obtain 



dN 




(49) 


At this order of the expansion, one can collect in (l49l l the zero-point contribution of the modes (the term 1/2 in 
between parentheses) and the bare chemical potential po,o(po) to form the true chemical potential poipo) of the fluid at 
zero temperature, and one can identify co.o in the prefactor of the phonon number operator BqBq with the true sound 
velocity at zero temperature, cq = [poPg(po)/nj]^^^- One then obtains the (low temperature) phononic limit of relation 
(Ell, without any constraint on the interaction strength. 


20. One first checks that the density fluctuations on a given lattice site ai'e small in relative value, using At T = 0, {6p^)lp^ ^ ^ 

using ou and the necessary condition mcllh > 1 established in note 19. At 7 > 0, there is the additional thermal contribution « ^ which 
is <s: 1 for the same reasons. Second, one checks that the phase fluctuations between two neighboring lattice sites are small in absolute value. To 
this end, one notes from |46) that /^((gradd^Z) = {mcllh)^{6p^)lp^. To conclude, it remains to use mcKhp^^^) = 0(1), a property that holds in 
the whole BEC-BCS crossover, as well as the previous estimates of {6p^)lp^. 
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Brouillage thermique d’un gaz coherent de fermions 
Thermal blurring of a coherent Fermi gas 
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Resume 

On considere generalement que la fonction d’onde macroscopique decrivant un condensat de paires de fermions possede une 
phase parfaitement definie et immuable. En realite, il n’existe que des systemes de taille finie, prepares, qui plus est, a temperature 
non nulle; le condensat possede alors un temps de coherence fini, meme lorsque le systeme demeure isole tout au long de son 
evolution et que le nombre de particules N est fixe. La memoire de la phase initiale se perd a mesure que le condensat interagit 
avec les modes excites, qui agissent comme un environnement dephasant. Get effet fondamental, crucial pour les applications qui 
exploitent la coherence macroscopique du condensat de paires, reste tres peu etudie. Dans cet article, nous relions le temps de 
coherence a la dynamique de phase du condensat, et nous montrons par une approche microscopique que la derivee de I’operateur 
phase du condensat 6q par rapport au temps est proportionnelle a un operateur potentiel chimique que nous construisons, et qui 
inclut les deux branches d’excitations du gaz, tant par brisure des paires que par mise en mouvement de leur centre de masse. 
Pour une realisation donnee d’energie E, 6q evolue aux temps longs comme ou fimc(E) est le potentiel chimique 

microcanonique; les fluctuations de I’energie d’une realisation a I’autre conduisent alors a un brouillage balistique de la phase, 
et a une decroissance gaussienne de la fonction de coherence temporelle avec un temps caracteristique oc En revanche, en 
I’absence de fluctuations d’energie, le temps de coherence diverge lineairement en A a cause du mouvement diffusif de Oq. Enfin, 
nous proposons une methode permettant de mesurer avec un gaz d’atomes froids ce temps de coherence, que nous predisons etre 
de I’ordre de la dizaine de millisecondes pour un gaz de fermions prepare dans I’ensemble canonique a la limite unitaire. 
Mots-cles : gaz de fermions ; fluides quantiques ; coherence quantique ; atomes froids 


1. Position du probleme 


Les gaz coherents d’atomes froids dans des pieges immateriels non dissipatifs sont des exemples uniques de 
systemes quantiques macroscopiques isoles. Se pose alors le probleme de leur temps de coherence intrinseque. Cette 
question fondamentale presente egalement un interet pratique, pour toutes les applications qui mettent a profit la 
coherence macroscopique, par exemple en interferometrie ou en ingenierie quantique oil Ton cree des etats intriques 
non triviau^ar evolution coherente Qui . Des mesures de temps de coherence sont deja accessibles avec des gaz 
de bosons Les experiences sur les gaz de fermions, qui se sont pour F instant concentrees sur les proprietes 

thermodynamiques Hi , done sur les aspects traditionnels du probleme a N corps, commencent a s’interesser aux 
mesures de correlations et de coherence |@], ce qui ouvriraun nouveau domaine de recherche, y compris dans le regime 
d’interaction forte : Foptique quantique fermionique Ii3]. Cependant, il n’y a a cejour aucune prediction sur le temps 
de coherence d’un condensat de paires de fermions, sauf a temperature nulle 111 111 . Dans cet article nous presentons 
la premiere theorie microscopique permettant de combler ce vide theorique en toute generalite. Notre analyse vaut 
pour d’autres systemes physiques, tels que les jonctions de Josephson mesoscopiques, pourvu qu’on arrive a reduire 
suffisamment la decoherence due a F environnement. 

Pour un gaz de bosons condense, le temps de coherence est determine par Fetalement de la distribution de probabi- 
lite du changement de phase du condensat pendant t. A temperature nulle et en presence d’interactions, cet etalement, 
balistique, est du aux fluctuations du nombre de particules. Cet effet a ete observe par interference de deux condensats 
initialement coherents, dont le nombre de particules fluctue du fait du bruit de partition 0i . Contrairement au cas 
du laser, qui est un systeme quantique ouvert, et de fa 9 on inattendue, un etalement balistique de la phase subsiste dans 
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ces systemes isoles pour un nombre de particules fixe si le gaz est a temperature non nulle lll2Lll3n . les fluctuations de 
cette autre quantite conservee qu’est I’energie rempla^ant celles du nombre de particules. 

Dans le cas d’un gaz non polarise de fermions condenses par paires, I’etude du temps de coherence presuppose 
une definition claire de la phase du condensat, et une construction explicite de I’operateur Oq correspondant ll ill . A 
temperature non nulle elle requiert en outre de pouvoir inclure dans la vitesse de variation de la phase la contribution 
des deux branches d’excitation, celle fermionique (par brisure des paires) et celle bosonique (par excitation du centre 
de masse des paires). Pour la branche fermionique, la RPA d’Anderson Cl suffit. Pour la branche bosonique en 
revanche, nous avons besoin d’un equivalent pour les fermions de la methode de Bogolioubov permettant de construire 
les operateurs de creation bl. et d’annihilation ba de quasi-particules et d’exprimer d0o/df en termes de ces operateurs. 
Qui plus est, il faut prendre en compte les interactions entre les quasi-particules dans revolution des ba- Ceci est une 
consequence non evidente du fait que la fonction d’onde du condensat de paires depend du nombre total de particules 
N, meme dans un systeme spatialement homogene, et va au-dela du programme de la RPA. 


2. Decroissance de la fonction de coherence 


Aux temperatures inferieures a la temperature critique, la fonction de coherence temporelle du champ de paires 
i^j(r)i^l(r'), ou i^o-(r) est I’operateur champ fermionique pour la composante de spin cr, est dominee aux temps longs 
par celle du condensat: 


giit) = {a'it)ao(Q)) 


( 1 ) 


ouao - f d^rd^r'(pQ(r, r')i^i(r)i^f(r') est la composante du champ de paires sur la fonction d’onde du condensat d- 
Comme le systeme a I’equilibre est dans un melange statistique d’etats propres a N corps avec des poids fl^, 
nous sommes ramenes a I’etude de la fonction de coherence g'l(t) dans I’etat propre d’energie £,) et de nombre de 
particules Na- Pour tirer parti des faibles fluctuations relatives du nombre de paires condensees pour un grand systeme, 
on utilise la decomposition d en les operateurs hermitiens phase et module 


ao 




( 2 ) 


et on approxime No, operateur nombre de paires condensees, par sa valeur moyenne No dans I’etat d’equilibre du 
systeme, pour obtenir 

gf(t) - (3) 

oil r operateur W, difference entre le transforme de H par e‘®" et H, 


fp = -H = -i[0o, H] - ^[00, [00, //]] + ... 


1 , 


(4) 


est approximativement N fois plus petit que H. En effet, comme ao, change le nombre total de particules d’un 
O(N^). Alors que H est une observable extensive, W est un operateur intensif et le double commutateur dans (lU est 
d’ordre 1/A. II apparait formellement dans I’equation (O un operateur d’evolution correspondant au hamiltonien H 
perturbe par W, restreint a I’etat propre de H. La fonction g */No est done, a un facteur de phase pres, proportion- 
nelle a I’amplitude de probabilite que le systeme initialement prepare dans y soit encore au bout d’un temps t en 
presence de W. Une approche naturelle pour obtenir une approximation non perturbative de cette amplitude est celle 

des fonctions de Green ou de I’operateur resolvante G(z) - (zl + lU)) du hamiltonien perturbe. Par la methode 
des projecteurs (voir la sous-section §III.B.2 de la reference OlSll l. on definit un hamiltonien effectif non hermitien 

Hes(z) gouvernant revolution restreinte a \t/fA), done tel que (iAi|Gfe)|iAi) “ “ (lA.il^eftWIiAi)) ■ Ceci conduit a 


siit) 


No 


r ^ 
Jc 2ot 


-Kz-Ei)tlh 


c2inz-{H^siz))A 


avec {H,s{z))a = + <W>i + {WQa 


Qa 


zQa-Qa(H+W)Qa 


QaW}a (5) 


Le domaine d’integration C dans le plan complexe est toute droite z = x + ir], rj > 0 fixe, decrite dans le sens des reels 
X decroissant de -i-oo a -oo, et Ton a introduit la notation {A)a = {il/A\A\^A) et I’operateur Qa - 1 — \^a){'4>a\, projecteur 
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orthogonal a En se limitant dans (i/eff(2))/i a I’ordre deux en W, done en omettant QaWQa au denominateur, et 
en negligeant sa dependance en z (approximation du pole), (Heff(z))i a; {H^fi(EA + iO^))i, on obtientQ 0 

gf(f) - avec tKSA-m)^{WQA - ^ (6) 

(Ea + i0+)2i - QaHQa 

On notera que 2ja n’est autre que le taux de depart de I’etat |i/r^) induit par W tel que le predit la regie d’or de Fermi. 

Le terme dominant sous I’exponentielle dans (|6]l est (lE),i, d’ordre comme W. Alin de I’interpreter, on constate 
d’abord que, d’apres le developpement (|4]i, en point de vue de Heisenberg, 


A ce stade il pourrait sembler difficile d’obtenir une expression explicite de d0o/df qui permettrait d’aller au-dela 
d’une description purement forme lie de la dynamique de phase. Heureusement il n’en est rien car, comme nous le 
montrerons dans la section suivante, la moyenne temporelle a gros grains de d0o/df dans un gaz faiblement excite est 
proportionnelle a un operateur potentiel chimique, quantite thermodynamique par nature : 


h d0o 
2"d7 


s=F,B O' 




( 8 ) 


La somme au second membre porte sur les deux branches d’excitation, celle des quasi-particules fermioniques (a 
inclut un indice de spin et un indice orbital, a = k, cr pour un systeme spatialement homogene) qui presente une 
bande interdite, et celle des quasi-particules bosoniques (a est uniquement orbital, a - q dans I’exemple precedent) 
dont le depart, a la limite thermodynamique et pour un systeme homogene, est phononique. Lorsque nous disons que 
le gaz doit etre faiblement excite, nous voulons dire que la depletion thermique du condensat doit rester petite. Ceci 
suppose en particulier que le nombre moyen de quasi-particules soit une faible fraction du nombre total de particules. 
La moyenne temporelle a gros grains est prise sur un temps long devant I’inverse des pulsations propres de ces 
quasi-particules mais court devant leur temps de collision, e’est-a-dire le temps typique de variation des operateurs 


1. Il y a den'iere 1’approximation du pole I’hypothese que le continuum auquel \{}/x) est couple est large en frequence (H . Aussi la loi exponen- 
tielle dans ID ne vaut-elle qu’aux temps longs devant I’inverse de cette largeur, e’est-a-dire ici devant le temps de correlation Tc des quasi-particules 
introduit plus bas. 

2. Pour faire apparaitre une ligne de coupure dans la resolvante puis un pole dans son prolongement analytique, done pour justifier rigoureu- 

sement I’emergence d’une partie imaginaire yx non nulle et donner un sens precis a I’equation il faut passer a la limite thermodynamique 
{N —^ +00 a densite fixee). Mais e’est ici inhabituellement subtil : comme le montre une generalisation du raisonnement autour de I’equation 
(To), la fonction deplacement Ka(z) = (WQa ^ a./i QaW)a, avec Imz > 0, tend vers zero en l/N, tout comme 6;i - iy^ dans (6). 

QAZ+QAiF^-H)Q^ 

L’astuce consiste alors a se placer a une echelle de temps d’ordre N, en posant t = Nr et en gai'dant r fixe a la limite thennodynamique. On 
elfectue le changement de variable z = + (W)/( + z' /N dans I’integrale de 0, on choisit t] = j]' jN {tj' > 0 fixe) dans le domaine d’integration 

C = {z = X + IT], X = +00 -oo}, ce qui conduit au domaine d’integration C' = [z' = x' + irf',x' = +oo -oo} sur z!, et Ton neglige le terme 

(^)/i fiui apparait au denominateur de la fonction deplacement ainsi que Re^'/A^ (par coherence avec I’omission de Qx^Qa). Le facteur 1/A dans 
I’element d’integration dz = dz'/N permet de soitir un facteur 1/A au denominateur de I’integrande, et Ton obtient a la limite thermodynamique le 
resultat suivant, independant de rj' par analyticite de I’integrande dans le demi-plan Imz' > 0 : 


lim 

N—*+oo 


Wo 


r dz' 

Jc' 2m z' - Mix 


avec Mix = lini AA^i{i? 7 VA) 

A^^-l-OO 


Voila qui foumit, pour un systeme de taille assez grande mais finie, la definition A(^^ - iy^i) = Si Ton avait pris la limite thermodynamique a 
t fixe avant de faire tendre t vers I’infini, on serait tombe, en developpant I’integrande de 0 au premier ordre en la fonction deplacement, sur une 
autre definition : lim/v^+oo A[Aq - 1] ~ (-iOA((5,i - iy^) avec N{6x - iyx) - 1111 ^ 77 ^ 0 + limA7^+co AA,i(i;;). Pour montrer I’equivalence 

de ces deux definitions, injectons dans la fonction deplacement une relation de fenneture sur les etats propres de H, A,i(z) = 'Zjfi:^x /(z + 

Ex - Efj), et voyons si la granularite de la distribution des - Ex dans cette somme ponderee tend vers zero plus vite que t]' /N. Pour cela, voyons 
chaque etat propre comme la superposition coherente d’un etat de Fock central de phonons et d’un faible halo d’etats de Fock de phonons s’en 
deduisant par des processus Beliaev-Landau 1 phonon 2 phonons. Meme si W conserve les nombres de phonons, cf. il peut coupler I’etat de 
Fock central de au halo de 1^^). Au pire, un seul processus Beliaev-Landau est mis en jeu, auquel cas E^ - Ex = ±(f 5 .ki + ^B.k 2 ~ ^ 5 .ki-i-k 2)5 
ki et k 2 sont les vecteurs d’onde des deux phonons emis ou absorbes ; comme la paire {ki,k 2 } prend oc valeurs differentes, ou V est le volume 
du systeme, E^ - Ex varie alors par pas oc 1/F“, ce qui est bien t]'/N. 
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nombres de quasi-particules hs^a, ce qui suppose que ces quasi-particules sont dans le regime faiblement collisionnel. 
Enfin fioiN) est le potentiel chimique a temperature nulle du gaz non polarise a N particules, autrement dit la derivee 
de I’energiedu fondamental par rapport a Nous interpretons le deuxieme terme du second membre de ([Sll comme un 
« operateur potentiel chimique » au sens ou sa moyenne quantique est la derivee adiabatique, c’est-a-dire a populations 
des quasi-particules fixees, de I’energie des quasi-particules Yjs=f.b Ha ^s,a{ns,a) par rapport a N. L’equation (O 
re lie la derivee de la phase au potentiel chimique au niveau des operateurs quantiques d’une theorie microscopique 
multimode. De ce point de vue, elle va au-dela de la seconde relation de Josephson etablie pour la phase du parametre 
d’ordre d’un supraconducteur (c/. la section §3.4 de la reference lllhll '). 

En prenant la moyenne de I’equation (|8ll dans \^,\) et en utilisant I’hypothese de thermicite des etats propres 
ifnl pour identifier moyenne dans un etat propre et moyenne d’ensemble microcanonique, on fait done apparaitre le 
potentiel chimique microcanonique/Tmc a I’energie et au nombre de particules Na '■ 


, )/l “ 2/Jn: 

df 


AEa,N,x) 


(9) 


Nous omettons la moyenne tempore lie a gros grains, ici superflue puisque la moyenne quantique est prise dans un etat 
stationnaire du systeme. 

Le terme suivant sous I’exponentielle dans (|6]l est d’ordre N * done sous-dominant. Pour le voir, on le relie a la 
fonction de correlation temporelle de d0o/df dans I’etat A une contribution d’ordre 1/A^ pres : 


r.i + = 



d6o(0 d0o(O)\ 
df df 



( 10 ) 


qui redonne bien (|6]l apres injection d’une relation de fermeture sur les etats propres de H. L’integrande a comme 
valeur initiale Vari(d0o/df) = (9(1/A) (ceci s’obtient en additionnant les variances des nombres de quasi-particules, 
ce qui correspond a un gaz parfait de quasi-particules dans 1’ensemble canonique et surestime done la variance micro¬ 
canonique) et decroit sur le temps typique de correlation Tc des hs^a, c’est-a-dire le temps typique de collision entre 
les quasi-particules. Ceci permet d’estimer |y,i -i- i(5i| a; TcV ar(d0o/df) = 0(1 IN). Le deplacement 6^ est done du 
meme ordre en N que le terme sous-dominant [0o, [^o, ^]] dans W ; tons deux sont N fois plus faibles que (|3) et on les 
neglige pour un grand systeme. En revanche ja est le seul a conduire a une decroissance exponentielle de la fonction 
de coherence microcanonique gj(f) et doit etre conserve. L’equation (fTOl i en fournit aussi une interpretation physique, 
si Ton se souvient que, dans la theorie du mouvement brownien, I’integrale de la fonction de correlation de la vitesse 
donne le coefficient de diffusion de la position : ja n’est autre que le coefficient de diffusion de la phase du condensat 
de paires lorsque le systeme est prepare dans I’ensemble microcanonique correspondant a ((Pa), et sera desormais note 
7 i = D(Ea,Na). Nous retenons flnalement : 

gflt) ^ 

L’equation (fTOl i fournit meme une methode de calcul de D : si Ton dispose d’equations cinetiques pour les nombres de 

quasi-particules apparaissant dans d0n/df, voir I’equation (|8ll, on pent calculer leur fonction de correlation temporelle, 
comme il est fait dans la reference llffll . et en deduire que 

D(E,N)^-(pX)-M-^C^^pA (12) 

La matrice M, de coefficients M^a.s'a', est la matrice des equations cinetiques linearisees = M6n donnant revo¬ 
lution des fluctuations diis^a des nombres de quasi-particules, regroupees en un vecteur 6n, autour de leurs valeurs 
stationnaires hs,a- Pour deflnir les autres notations, on introduit comme dans Q la matrice Ccan des covariances des 
nombres de quasi-particules „ dans 1’ensemble canonique a N particules et d’energie moyenne E, le vecteur energie 
e de composantes ejQ. et son vecteur dual e oc Ccanf normalise defa 9 on que 1.0Alors le vecteur A = (llftfAejAN 


3. On a (Ccan).^,^'^' = ± ns,a) ris,a = ^I^bT) t 1] OU le signe superieur (inferieur) vaut pour la branche d’excitation 

bosonique (fermionique). ^ et ^ sont des vecteurs propres a gauche et a droite de M avec la valeur propre nulle, comme le montrent la conservation 
de I’energie et I’ecriture plus parlante ?oc dn/dr, T etant la temperature de I’ensemble canonique et n le vecteur de composantes hs^a S. 
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regroupe les coefficients des hs,a dans I’expression ([ 8 ]) de -dOo/dt , la matrice P est le projecteur non orthogonal^ tel 
que P6n - 6n- e{e ■ 6n) pour tout 6n, et la matrice des covariances microcanonique vaut Cmc = P^ CcmP 

II nous reste maintenant a prendre la moyenne statistique de 1’equation (HUl sur les etats constituant I’etat 
melange du systeme. Aux grands N, on suppose que I’ensemble generalise donne par les Ef^ conduit a des fluctuations 
de I’energie totale et du nombre de particules tres faibles en valeur relative autour de leur valeur moyenne E et de 
N (c’est le cas par exemple des ensembles canonique et grand canonique). On en prend done une approximation 
gaussienne et on linearise /ime autour de {E, N). A cet ordre, on pent remplacer D, deja sous-dominant d’un facteur N, 
par D{E,N). II n’y a plus qu’a calculer I’integrale gaussienne du facteur de phase exp{2i[(diJ.„iciE,N)/dE){E - E) + 
{dfimciE, N)/dN){N - N)]t/h} pour obtenir le resultat central de ce travail: 

gl{t) (I 3 ) 

En presence de fluctuations de I’energie ou du nombre de particules, le brouillage thermique aux temps longs prend 
la forme d’une decroissance gaussienne de la fonction de coherence avec un temps caracteristique 

= Var|A^(£,A) + £^(£,A)j (14) 

qui diverge en pour des fluctuations normales. Le coefficient de diffusion de la phase D donne, quant a lui, un 
temps de decroissance caracteristique qui diverge lineairement en N ; c’est done bien une contribution sous-dominante 
aux temps longs, sauf dans I’ensemble microcanonique ou cette diffusion de phase intrinseque pent etre observee 
directement. 


3. Derivation microscopique de I’equation sur I’operateur phase 


Comme promis, nous donnons maintenant a notre connaissance la premiere demonstration microscopique de 
I’equation ([ 8 ]), reliant la vitesse de variation de I’operateur phase d’un condensat de paires a ce que Ton peut qualifier 
d’operateurpotentiel chimique. 

La contribution de la branche d’excitation fermionique a d0o/df peut s’obtenir par des equations du mouvement 


linearisees pour de faibles fluctuations des operateurs de paires et autour de I’etat de champ moyen, 

ce qui constitue la Random Phase Approximation (RPA) d’Anderson il41] . En effet, en calculant la moyenne tempo- 
relle de dOo/dt grace a 1 ’equation ( 120 ) de la reference 111 If et en recrivant I’equation ( 86 ) de cette reference en termes 
des nombres d’occupation des quasi-particules fermioniques,]^ on obtient 


h d0o 

2"dr 


RPA^ dN 


.A)+ X 

Q'=k,c 


defo 


dN 


- f^F,a 


(15) 


ou efXa- est le spectre d’excitation BCS, ici pourun systeme spatialement homogene, et N le nombre moyen BCS de 
particules dans I’ensemble grand canonique de potentiel chimique /i .3 

Des difficultes fondamentales nous empechent d’obtenir la contribution de la branche de phonons a F equation (l 8 ]l 
dans le cadre de la RPA. 0 Nous avons alors choisi de traiter le probleme par une approche variationnelle avec comme 


4. P? = 0 et ft projette sur le sous-espace des fluctuations dn d’energie nulle, {6n | e' ■ dn = 0|, a I’interieur duquel I’inverse de M est defini. 

5. On remarque que +flF,-k,i oil le membre de gauche fait reference aux notations de I’article (H , et on se sert de son equation 

(74) pour reconnaitre d6/rk,o-/d/i. 

6. L’equation (B) est ecrite en reintroduisant le terme de phase trivial -2fi(N)tlh qui n’appai'ait pas dans la reference [H du fait de I’utilisation 
du hamiltonien grand canonique. 

7. Le resultat GD restreint a I’ordre lineaire en les operateurs de paires n'inclut pas la contribution de la branche de phonons. On pourrait 
esperer obtenir cette contribution en poussant le calcul par la RPA jusqu’a I’ordre quadratique comme suit : premierement, on calcule dOo/dt 
jusqu’a I’ordre quadratique. Deuxiemement, on forme I’ensemble des operateurs de la RPA qui transferent une quantite de mouvement hq, a savoir 

^k+qi^-ki’ ^k+qT^’^f ^^ihilc un fermion de vecteur d’onde k et de spin f, ensemble qui est clos au regard des 

equations lineaiisees. Pour chaque valeur de q la matrice Xq de ce systeme d’equations possede deux modes collectifs, dont les energies, opposees 
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Ansatz I’etat coherent de paires dependant du temps le plus general^ 


|iA) = A^(f)exp 


Z7®^r(r,r';r)(Aj:(r)^|(r') 


r.r' 


10 ), 


(16) 


Cet Ansatz, fonction des variables independantes r(r, r') et normalise a I’unite grace au facteur N{f), est ecrit ici 
apres discretisation de I’espace reel en un reseau cubique de pas I (que Ton fait tendre vers 0 a la fin des calculs), 
les operateurs champ obeissant alors a des relations d’anticommutation du type {i^o-(r), = bo-.o-'br.r'/l^- Lejeu 

de variables canoniquement conjuguees associe, 0(r, r'), ®*(r, r') est construitl^ dans la section §9.9b de la reference 
S. II ne faut pas confondre le champ <5 avec le champ de paires (iAx*At) [3 Quand les paires de I’etat (fTbl l sont au 
repos, r(r, r', t) depend uniquement de r - r' et la transformee de Fourier de <I> n’est autre que I’amplitude BCS Vk de 
la paire k t, -k i in , tandis que la transformee de Fourier du champ de paires s’ecrit -t/kFk = -Fk(l - |FkP)0^. 
Lorsque les paires sont en mouvement, nous n’avons plus d’interpretation physique simple mais la norme au carre de 
O demeure egale a la moitie du nombre moyen de particules dans I’etat |i/r): 

^ = ||(b||2 = /®^|<l)(r,r';0l^ (17) 

r,r' 


Le hamiltonien regissant revolution du champ d), a savoir 

-7/(0), <b*) = 


(18) 


pourrait etre explicite a I’aide du theoreme de Wick, mais nous aurons ici besoin seulement de son invariance par 
changement de phase globale 4>(r, r') —> e‘^<l>(r, r'), Vy e R (symetrie t/(l)), qui decoule de la conservation du 
nombre total de particules N par H. A temperature nulle, le champ <l>(r, r') a N fixe est fige, a une phase globale pres, 
dans le minimiseur <l)o(r, rO = (lV/2)O^0o(r, rO de 'H, oil 0o, reel et normalise a I’unite, depend de N meme dans le 
cas spatialement homogene, differe de la fonction d’onde du condensat de la meme maniere que O differe du 
champ de paires (voir la note 10). A temperature suffisamment faible, on pent developper”?/ en puissances des 

petites deviations de d) par rapport au cercle y i-» e‘r(I)o(r, r'), lieu des minima de "Tf a A fixe. On decompose done le 
champ en ses composantes colineaire et orthogonale a (po : 

0(r, r') = e‘®[n‘^Vo(r, r') + A(r, r')] (19) 


Tune de I’autre, sont lineaires en q dans la limite des grandes longueurs d’onde. On developpe alors les operateurs de la RPA sur I’ensemble des 
modes propres de Xq ; les amplitudes obtenues sur les modes collectifs, ^_q et ^q, annihilent et creent respectivement une excitation collective 
d’impulsion Thq. Troisiemement, on insere cette decomposition modale dans les operateurs de la RPA apparaissant dans dOo/dt pour obtenir une 
expression quadratique en les amplitudes modules, et done des termes en b^bq comme dans m- Malheureusement, cette approche par la RPA n’est 
pas fiable et nous devons y renoncer car les operateurs de la RPA, bien que lineairement independants, ne sont pas quadratiquement independants, 
comme on le voit en dissociant et reassociant les operateurs de paires grace aux relations d’anticommutation fermioniques, par exemple 

(«LqT«kT)(alk-qi“-ki) = («LqT“-k-qi)(«-kAkT). 

Ceci montre (/) qu’il n’y a pas une maniere unique d’exprimer dOo/dt comme une fonction quadratique des operateurs de la RPA, (n) que les 
amplitudes modules ne sont pas quadratiquement independantes, et (Hi) que le coefficient de b^bq n’est pas determine de maniere unique par cette 
approche fondee sur la RPA 

8. Nous utilisons ici par souci de simplicite un Ansatz de type BCS dependant du temps dans lequel le nombre de particules possede des fluc¬ 

tuations quantiques (a ne pas confondre avec les fluctuations thermiques de I’ensemble grand canonique). Au prix d’une difficulte supplementaire, 
nous aurions pu tout aussi bien choisir un Ansatz a nombre de particules fixe fT3l . ce qui n’aurait change les coefficients de la fonctionnelle energie 
que d’une correction relative si bien que le spectre obtenu a la limite thermodynamique eut ete le meme. 

9. Nos variables F et <!> correspondent aux variables z et yS de la reference fl^ dans le cas ou I’on remplace les operateurs et le vide de quasi- 

particules dans I’equation (9.132) de (H par ceux des particules. Si E et ^ sont les matrices de coefficients respectivement /^r(r, r') et r') 

alors ^ = -E(l + d’apres I’equation (9.144) de fT^ . Leurs equations d’evolution s’obtiennent par methode vaiiationnelle a partir du 

lagrangien L = ih [(j/fl(d/d?)|(/r) — c.c.] /2 - (^|//|</f). <1) a ete choisi astucieusement de sorte que L = Sr,r' r')<9r<I>(r, r') - c.c.] /2 - 

7F(<I), O*), ce qui conduit a dt^ = . Ceci montre que la vaiiable conjuguee de <1) dans le formalisme hamiltonien est pour le 

crochet de Poisson traditionnel, et done pour le crochet de Poisson {C>, O*} = (i^/^)“^, sachant qu’on doit avoir diO = {<I),7F}. 

10. Avec la meme notation matricielle que dans la note 9, nous avons (^-ft^j^) = ^(1 - d’apres I’equation (9.146) de la reference IT3l . 
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ou I’excursion de la phase 6 peut etre arbitrairement grande, contrairement a celle de A. Ce cadre permet de deve- 
lopper une theorie perturbative systematique en puissances du champ A (voir I’annexe A), qui est I’analogue pour les 
fermions de la theorie de Bogolioubov sans brisure de symetrie t/(l) pour les bosons 12(111 . Comme nous le verrons, 
la phase 0 reste proche de la phase Oq du condensat si A est assez petit. On ecrit done les equations du mouvement 
pour 6 et pour les champs de I’espace orthogonal A, A*. Dans la forme finale des equations, on fait systematiquement 
disparaitre les variables du condensat grace a la relation n - N/2 - ||A|p, consequence de (fT71) . et on se limite a I’ordre 
2en A, A*. 

La difhculte centrale de notre probleme provient de F existence d’un terme lineaire en A, A* dans d6/dt, conse¬ 
quence de la dependance de en le nombre de paires im . En r absence de ce terme lineaire, il suffirait de developper 
le champ A sur les modes propres de ses petites oscillations lineaires, deduits de F approximation quadratique du ha- 
miltonien “K a A fixe : 


/A(r,r';f)\ 


Yjbait) 


luair, r')\ 
r')/ 


+ Kit) 


A'a(Lr')'l 


( 20 ) 


ou la somme porte sur les modes propres d’energie positive e„, avec la normalisation ||Ma|p - ||Va|P = 1. A cet ordre, 
bait) = baiQ)e^"^‘‘‘^^■ II n’y aurait qu’a injecter la decomposition (l20l i dans dO/dt et a prendre une moyenne temporelle 
a gros grains pour tuer les termes de battement entre les modes; dO/dt contiendrait alors la combinaison lineaire 
attendue des nombres de quasi-particules bosoniques La realite est plus subtile : a cause de Finteraction 

entre les quasi-particules, ba n’est pas nul mais d’ordre deux en A et A*, si bien que le terme lineaire dans dO/dt 
contribue a la meme hauteur que les termes quadratiques. On doit done calculer ba , ce que Fon fait dans Fannexe 
A, en utilisant en particulier le caractere borne du champ A, consequence du respect de la symetrie U{1) dans le 
developpement (fTOT l, puis en appliquant le theoreme de Hellmann-Feynman, pour trouver : 


h dQ 
2 df 


tiQim + Y, + odiAih 

a 


( 21 ) 


Discutons brievement le spectre dans le cas d’un systeme spatialement homogene, a la limite continue I —» 0 pour 
une interaction de contact dans Fonde s de longueur de diffusion fixee entre fermions de spin oppose. Pour chaque va- 
leur du vecteur d’onde total q, il se compose d’au plus une energie discrete et d’un continuum parametre par deux 
vecteurs d’onde (ki, t; k 2 , i) eF.k,,T + fF.k 2 ,i de somme fixee (ki H-k? = q), ouef t.o- est la relation de dispersion BCS. 
La branche q i-> q coincide bien avec celle de la reference 12111 . comme nous Favons verifie. De depart phononique, 
elle correspond aux excitations elementaires bosoniques du gaz de fermions dont nous cherchions la contribution a 
la dynamique de phase. Elle doit done etre gardee dans (ISTTi . En revanche, le continuum correspond a Fexcitation de 
deux quasi-particules fermioniques; celles-ci ne peuvent en effet etre excitees que par paires a partir du fondamental 
car le hamiltonien H contient un nombre pair de facteurs \j/ qui detruisent ou creent chacun une quasi-particule. 
Ces biexcitations ne sont pas physiquement independantes,[3 et font double emploi avec la contribution de la RPA a 
d0o/df. Il ne faut done pas les inclure dans (1211 1. 

Deux dernieres remarques nous separent de (O. (/) Les champs ) et <1) sont distincts, aussi les phases 0o et 6 
ne coincident-elles pas exactement. Leur valeur moyenne temporelle a gros grains ne different cependant que par un 
terme d’ordre ||A|p, borne, negligeable aux temps longs, ne contribuant pas au brouillage de phase du condensat de 
paires. En i) La phase de notre approche variationnelle est une variable classique, alors que Pq dans (|8]l est un operateur 
quantique. Ce hiatus est comble par application au resultat (l2T1 i de la procedure de quantification du chapitre 11 de la 
reference Cl , puisque les ba de la branche bosonique y sont remplaces par des operateurs bosoniquesl^fio,, [ba, K] - 


11. Exciter O' = (ki, t;k 2 , i) et ff' = (k'j, t; kj, i) revient a exciter O'" = (k'j, t; k 2 , i) et cr'" =(ki,t;kLi). 

12. En exprimant dans (^o) la valeur moyenne du champ de paires en fonction de O, on se rend compte pour A faible que {ao) = V^(l + 
0(||A||)) et done, etant donne que A = 0(||A|p), on a 0o =6 + 0(||A|p). 

13. Plus precisement, ces operateurs ne sont bosoniques que pour une faible densite d’excitations. Pour un systeme homogene spatialement et 

dans un formalisme qui brise la symetrie t/(l) {6 = 0), nous obtenons a partir des equations (11.81c) de fT^ . generalisees au cas de fermions 
apparies sont remplaces par les operateurs BCS ^ko- d’annihilation de quasi-particules fermioniques) et restreintes aux images bosoniques 

faiblement excitees (BB^ negligeable) : = Zk ^k^-ki^k-fqt + ^k^I-qt^-ki' coefficients reels et sont des combinaisons lineaires des 

coefficients correspondants Wq et Vq dans I’equation et heritent de leur condition de normalisation Zk(-^k^^ ” 

avec e = 0, 5q = [fcq.feq] - 1 a une valeur moyenne et une variance 0(||A|p/Al) puisque &ko-|iA) = 0(||A||/Al'^^) si ||A||/A'^^ -*0. 
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Figure 1: Quelques considerations experimentales. (a) Configuration de piegeage proposee pour mesurer gi(0 par interferometrie de Ramsey : le 
condensat de paires de fermions est confine dans le piege principal (a fond plat sur la figure); on transfere au plus deux atomes en moyenne (sous 
forme d’un dimere) dans le piege secondaire (tres etroit) par effet tunnel resonnant d’amplitude ajustable pai' une bairiere de hauteur variable ; 
on cree ainsi une reference de phase avec laquelle on fait interferer le condensat de paires apres un temps d’evolution t. (b) Temps de brouillage 
thermique pour un gaz coherent de fermions a la limite unitaire dans 1’ ensemble canonique, en fonction de la temperature T en unites de la 
temperature de Fermi Tf = eflkg. Disques : a paitir de I’equation d’etat mesuree dans la reference a. Ligne tiretee : formule (22) qui utilise 
requation d’etat approchee decrite dans le texte. 


1. Nous pensons que la validite de I’equation (O, reliant d0o/df a I’operateur potentiel chimique, et done du resultat 
(foi l, depasse celle de notre derivation microscopique variationnelle et inclut le regime d’interaction forte pourvu que 
les quasi-particules aient une duree de vie beaucoup plus longue que I’inverse de leur frequence propre, ce qui est le 
cas a temperature sufSsamment basse. Dans le cas limite ou Ton pent negliger la branche d’excitation fermionique 
et se limiter a la partie phononique de la branche bosonique, e’est bien ce que predit la variante irrotationnelle de la 
theorie de I’hydrodynamique quantique de Landau et Khalatnikov (voir I’annexe B). 

4. En vue d’une etude experimentale 

Esquissons maintenant un protocole de mise en evidence experimentale du brouillage thermique d’un condensat 
de paires. L’idee centrale est, dans les phases de preparation et de mesure, de bosoniser de maniere adiabatique 
reversible les paires de Cooper atomiques sous forme de dimeres fortement lies et faiblement interagissants, en 
ajustant la longueur de diffusion a une faible valeur positive grace a une resonance de Feshbach magnetique. Ceci 
permet (/) de produire un echantillon a nombre faiblement fluctuant de dimeres, en faisant fondre une phase de Mott 
preparee par realisation du modele de Hubbard bosonique 1^ . (ii) de creer un couplage tunnel entre le piege principal 
renfermant les N particules et un piege secondaire tres etroit, a travers une barriere de hauteur ajustable yO (FiglT^), 
(Hi) de detecter par fluorescence un dimere unique lEsIl dans ce piege secondaire. Pour mesurer la fonction gi(f), nous 
adaptons aux paires de fermions la methode interferometrique de Ramsey de la reference lE^ . ou sont appliquees deux 
impulsions de Rabi a un intervalle de temps t. On prepare initialement les paires bosonisees dans le piege principal. 
La premiere impulsion est d’angle e suffisamment faible pour transferer en moyenne moins d’un dimere dans le piege 
secondaire et eviter ainsi le bruit de partition qui masquerait le brouillage thermique. On ajuste ensuite les interactions 
a la valeur a laquelle on souhaite etudier la dynamique de phase du condensat de paires, et on laisse evoluer pendant 
le temps t. On rebosonise et on applique une seconde impulsion d’angle e avant de mesurer le nombre «sec de dimeres 
dans le piege secondaire. Alors la valeur moyenne de nsec sur les realisations presente, en fonction du temps f, des 
oscillations de demi-pulsation la difference (sur h) des potentiels chimiques entre les deux zones de piegeage, et de 
contraste |gi(f)/.gi(0)|. 

5. Quelques resultats explicites sur le gaz unitaire 

Pour commencer, donnons une estimation du temps de brouillage pour un gaz unitaire prepare dans F ensemble 
canonique, done avec des fluctuations d’energie de variance VarF = kBT^drE. En utilisant I’equation d’etat du gaz 
unitaire non polarise mesuree dans la reference JH] nous trouvons, pour un systeme spatialement homogene (dans un 
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piege a fond plat 0271] ') les temps de brouillage thermiques fbr representes par des disques sur la figure[T} 5 . Ainsi, a une 
temperature T — 0,12 Tf 0,7 Tc ou Tc est la temperature de condensation des paires, on trouve fbr ~ er, ce 

qui correspond a 20 millisecondes pour une temperature de Fermi typique Tf — eflks - 1//K et un nombre d’atomes 
typique N - 10^. On pent aussi comme dans la reference estimer I’equation d’etat du gaz unitaire a partir d’un 
modele de gaz parfait de quasi-particules avec des relations de dispersion simples. Pour la branche bosonique, on prend 


1 2^ ffi.q = ficq avec c la vitesse du son aT - 0, mc^ - et ^ le parametre de Bertsch. Pour la branche fermionique, 
on prend 12^ ff.k.o- = A+(^^-eo)^/(2/o), oil A est le gap et ou eo et/o donnent la position du minimum et la courbure 
de la relation de dispersion. En gardant pour chaque branche sa contribution a I’energie volumique moyenne d’ordre 
dominant a basse temperature EfjV - 




et EbIV ^ ^{keTfUncf comme dans 0281], et en 


utilisant les valeurs experimentales ^ = 0,376, A = Q,44ef, eo - 0,85 €f et la valeur theorique ll^ /o = 0,846ef, 
nous trouvons ; 


(tbr^p)^ 
\9/2 


0 

0 ^ 


(1 + Irf 
(1 + f) 


( 22 ) 


oil & - T/Tp et oil r= dpEpIdpEB - e 0.44/0 mesure Fimportance relative des deux branches d’excitation. 

Cette formule, representee en tirete sur la figure[T] 3 , donne un equivalent exact de fbi quand 0 ^ O.Pl 

Aux temps finisEl t = <s; fbr, la contribution de D a gi (f) dans Fequation (fTSl l est a priori comparable a celle 

de puisque les deux coefficients se comportent comme 1/N dans Fensemble canonique. D pent etre calcule a tres 
basse temperature en ne gardant que la partie phononique du spectre. Dans I’etat de Fart, il est predit par differentes 
approches que la branche bosonique est convexe aux faibles q dans un gaz unitaire j 31, 32, 3, 




■ 0{q^) avec 7 > 0 


(23) 


si bien que les phonons interagissent par des processus Landau-Beliaev 2 phonons <-> 1 phonon comme dans un gaz 
de bosons en interaction faible.0 Nous pouvons alors reutiliser, en Fadaptant,0 la reduction a basse temperature de 
Fexpression generale (fTSl i effectuee dans la reference 11811 . D’autre part, nous trouvons, par resolution numerique des 
equations de la RPA de la reference 12111 , que y ^ 0,1, en accord avec il33Ll34ll . si bien que 


hND 


C0^ 




(24) 


6/r 0^0 

Ceci acheve notre description du brouillage thermique d’un gaz de fermions a la limite unitaire et a basse temperature. 


6. Conclusion 

Nous avons presente la premiere theorie microscopique du brouillage thermique de la phase d’un condensat de 
paires de fermions O, et mis en evidence un brouillage balistique et une diffusion de phase sous-dominante. Le 


14. Le bon accord avec les donnees experimentales est a relativiser car si Ton traite les deux branches du modele a tous les ordres en kgT, on 

obtient un decalage vers le haut de a peu pres constant et egal a 5 sur I’intervalle de temperature de la figure Et. 

15. Les temps que Ton considere sont neanmoins contraints par I’inegalite t > Tc (voir la note 1) : pour un gaz unitaire Tc-( 0, 149/0)^ quand 
0 —» 0 , d’apres la noteflT] 

16. Dans le cas concave, les processus dominants sont les diffusions 2 phonons 2 phonons comme dans la reference (22l 

17. L’astuce universalisante a basse temperature est d’utiliser, quel que soit le systeme physique, mc/h comme unite de vecteur d’onde et mc^ 
comme unite d’e nerg ie. Partons alors du coefficient de diffusion de la phase et du temps de con'elation des quasi-particules pour des bosons en 
interaction faible fl^ : hND/{mc^) ~ ciihgTjmc^)^ et mc^rdh ~ C3p[h/(^/2mc)f(mc^/kBT)^, ou c\ ^ 0,3036et C 3 ^ 0,05472, etrecensons tous 
les facteurs correctifs faisant passer du gaz de bosons en interaction faible au gaz unitaire de fermions, compte tenu du fait que D est une function 
quadratique de PA et lineaire de M~^ , voir I’equation {\2\ : (f) dans le cas d’un condensat de paires, il y a un facteur deux dans les coefficients A de 

-diSoldt , done un facteur quatre sur D, (ii) fequation d’etat passe de poip) oc papoip) oc p^^^, ou p est la densite, done (A^/c)dc/dA^ dans A passe de 
1/2 a 1/3, d’ou un facteur 2/3 sur A et un facteur 4/9 sur D, {Hi) PA, nul pour une branche d’excitation bosonique lineaire en q puisqu’alors A oc 
est proportionnel a basse temperature au parametre de courbure sans dimension 7 , qui valait un dans (H, d’ou un facteur 7 ^ dans D, (fv) d’apres 
I’hydrodynamique quantique (23l . les amplitudes adimensionnees de couplage Beliaev-Landau entre les quasi-particules, a des valeurs fixees des 
vecteurs d’onde adimensionnes, dependent du systeme seulement via un facteur global 1 + A, ou A = (p)/[3/./q(p)] passe de 0 a A,, = -1/9, 

d’ou un facteur (1 + Au)~'^ sur la matrice done sur D et sur le temps de correlation Tc induit par les collisions entre les quasi-particules. Alors, 
pour le gaz unitaire, hDNl{mc^) ~ ci[( 47 / 3 )^/(l + Au)^]{kBT/mc^)‘^ et mc^rdh ~ C3(l + A„)“^p[/i/( V2mc)]^(mc“//:sr)^. 
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temps de brouillage depend de la variance de I’energie totale du gaz, ainsi que de la derivee du potentiel chimique 
microcanonique du gaz par rapport a I’energie. Le coefficient de diffusion de la phase pent etre deduit d’equations 
cinetiques sur les nombres de quasi-particules. Pour etablir ces resultats, nous avons utilise de fa 9 on centrale le fait 
que la vitesse de variation de la phase du condensat de paires est reliee a I’operateur potentiel chimique du gaz par 
r equation ([8]), dont nous avons donne une derivation microscopique incluant les deux branches d’excitation du gaz, 
bosonique et fermionique. Nous avons enfin propose un protocole experimental de mesure de ce temps, que nous 
avons estime etre de I’ordre de la dizaine de millisecondes pour un gaz coherent de fermions prepare a la limite 
unitaire dans I’ensemble canonique. 
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Annexe A. Complements sur le calcul variationnel avec un etat coherent de paires 

Dans cette annexe nous demontrons F equation (l2Tli portant sur la moyenne temporelle de 9 grace a un modele 
microscopique fonde sur F Ansatz (fTbl l. etat coherent de paires en mouvement. II est commode dans un premier temps 
de traiter la variable reelle ii et le champ complexe A dans ( fT9b comme des variables independantes, c’est-a-dire de 
ne pas fixer ||<1)||. Pour inclure les interactions entre les quasi-particules, il faut aller jusqu’a Fordre trois en A et A* 
comme suit: 

3 

w, CD*) = To[n, <Po(N)] + ^ Tj[n, 0o(A)](A, A*) + D(||A|h (25) 

7=1 

Le tenseur Tj est de rang j, c’est-a-dire que Tj{K, A*) est d’ordre exactement j en A et A*; il pent s’exprimer en 
termes de la differentielle j-ieme de "TY prise en (0,0*) = restreinte au sous-espace orthogonal a 

(00,0) et (0,0o), ou Fon a pris 0o reelle, et il ne depend pas de la phase 9 du fait de la symetrie t/(l). A nombre total 
de particules fixe, Fenergie ne varie pas au premier ordre autour du minimiseur si bien que 71 [A/2,0o(A)] = 0. Cette 
innocente equation n’est autre que la celebre equation du gap lorsque le systeme est spatialement homogene. De plus, 
on verifie que (9„7o[A/2,0o(A)] = 2po{N) ou pq{N) = d£’o(A)/dA est le potentiel chimique du gaz a temperature 
nulle, Eq(N) — 7o[A/2, 4>q(N)] etant Fenergie de Fetat fondamental. 

Les variables phase et module au carre de Famplitude du champ O dans le mode 0o sont canoniquement conju- 
guees, si bien que -tid9fdt = 5„7Y(0, O*). Une fois effectuee cette derivee dans (l25l l a A et A* fixes, on pent fixer la 
norme de O a (A/2)' (c’ est-a-dire le nombre total de particules a A) et eliminer n grace a la relation ii = ||0|p-||A|p; 

le champ A reste alors la seule variable dynamique du probleme. L’expression est utile seulement jusqu’a Fordre deux 
en A, A* : 


“ " 5«^o[A/2, 00(A)] - \\AfdlTo[N/2, 0o(A)] + ^ d„Tj[N/2, 0o(A)](A, A*) + C>(||A||3) (26) 

7=1 

Le hamiltonien regissant Fevolution de A a nombre de particules fixe s’obtient en rempla 9 ant ii par A/2 - ||A|p 
dans (l25l l et en developpant jusqu’a Fordre trois en A, A* ; 

ElNiKA*) = £o(A) + f2[A](A, A*) + 7-3[A/2,0 o(A)](A, A*) - ||A||25„7-i[A/2,0o(A)](A, A*) + C>(||A||4) (27) 


ou la forme quadratique T 2 [A](A, A*) est obtenue en soustrayant 2/7o(A)||A|p a T ^[A/2,0o(A)](A, A*). Pour calcu- 
ler la moyenne temporelle a gros grains de 5„7“i[A/2,0o(A)](A, A*), c’est-a-dire du terme lineaire en A, A* posant 
probleme dans dP/df (voir le corps de Particle), il est astucieux d’ecrire alors directement la derivee temporelle de la 
partie imaginaire du champ A projetee sur la fonction (A/2)'^^d0o/dA, 



y/Ay/" d0o(r,rO 
A\2 / dA 


(A(r, r') - A*(r, r')) 


(28) 
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Comme ihBtK - I A*), il vient 


dF 

- 2h— = D ■ -KwCA, A*) = D ■ T2[N]{K, A*) + 2) ■ A*) + (9(||A|h 

df 


oil 'K™*’ est la partie de 'Hn d’ordre trois en A, A*. On a introduit I’operateur differentiel 

,1/2 


® “ Zj i 2 j dA '^^•(r.r')) 

r.r' 


( 29 ) 


(30) 


Nous allons maintenant mettre a profit en cascade deux relations verifiees pour tout A orthogonal a cpoiN ); 

22)-f2[A](A,A*) = -5„7-i[A/2,0o(A^)](A,A*) (31) 

2D-9{f(A,A*) = 2^f2[A](A,A*)-(i„T2[A/2,0o(A^)](A,A*) + ||A||2527-o[A/2,,^o(A^)] (32) 

On etablit ces relations en considerant formellement un champ O' a A + 6N particules et en obtenant de deux fa 9 ons 
differentes Fapproximation quadratique en A' et A'* de 'HN+sNiA', A'*) oil A' est a une phase pres comme dans (fT^ 
(ecrite pour A + 5N particules) la composante de O' orthogonale a 4>(i{N + 5N). D’une part, on remplace simplement 
A par A + 6N et A par A' dans (l27l i. puis on developpe au premier ordre en 6N, ce qui fait apparaitre dT 2 /dA[A]. 
D’autre part, on applique a '72(0', O'*) le developpement (l25T l autour de Oo(A) sachant qu’au premier ordre en 6N, la 
composante de O' orthogonale a (poiN) contient, en plus de celle A'j^ de A', une contribution venant de dipoldN : 

A(r, r') = (f ^(r, O + Al(r, r') + 0(6N\ ||A'||^ dA||A'||2) (33) 

C’est ce decalage inhnitesimal proportionnel a 6N dans la direction de d^o/dA qui fait apparaitre Foperateur 2). 
Les equations (EB et (132b resultent de Fidentihcation aux ordres un et deux en A' et A'* des deux expressions de 
^n+6n{A', A'*) obtenues. 

Le plus gros du travail a ete accompli. II reste a combiner les equations (|2^ . (l29l l. (l3TT i. et (l32l i. ou plus precisement 
leur moyenne temporelle a gros grains (sur une echelle de temps longue devant la periode d’oscillation des amplitudes 
modules ba mais courte devant le temps d’evolution des nombres de quasi-particules \ba\^) pour obtenir 


=|io(A) + ^[A](A,A*) +0(||A|p) (34) 

2 df dA 

ou nous avons utilise de fa 9 on cruciale la nullite de dA/df (et done de dF/df) pour le champ A, dont F excursion est 
bornee. La forme quadratique TtiN] est representee par la matrice rjXXN] avec 'f - soit, en notation par 

blocs et avec le produit scalaire hermitien (,) dont decoule la norme || || : 

■f2[A](A, A*) = i |(a a*) , l^J:m (35) 

II reste a injecter la decomposition module (l20l i dans Fequation dLSl l prealablement derivee par rapport a A, les 
et (Mq., Vq,), —ea et (v*, «*) etant les valeurs propres et vecteurs propres de La valeur moyenne a gros grains 

se charge d’eliminer les termes croises et le theoreme de Hellmann-Feynman assure que[3 {{ua,-Va), “ 

deo/dN, et permet ainsi de conclure sur Fequation (I 2 TI 1 . 


18. (Ua, -Va) est le vecteur dual de {Ua, Vq.) car = 7]-C- 
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Annexe B. Hydrodynamique quantique irrotationnelle 

Dans cette annexe nous donnons une demonstration supplementaire de I’equation Q a partir de Fhydrodynamique 
quantique. Bien que cette demonstration ne soit pas microscopique, a F inverse de celle que nous avons donnee dans le 
corps de Farticle et dans Fannexe A, et bien qu’elle neglige les degres de liberte fermioniques internes et se contente 
d’une description aux grandes echelles spatiales des paires de fermions, par un champ bosonique done, elle presente 
Favantage de reposer uniquement sur Fequation d’etat et d’etre ainsi valide dans tons les regimes d’interaction. 

En 1949, pour calculer la viscosite de F helium superfluide a basse temperature, Landau et Khalatnikov ont deve- 
loppe la theorie de F hydrodynamique quantique J^, qui permet de determiner, de maniere exacte a Fordre dominant 
en r, les effets d’une temperature non nulle sur les observables d’un fluide quantique, du moins sur celles dependant 
seulement des faibles echelles d’energie et des grandes echelles de longueur. De fa 9 on remarquable, le seul ingredient 
specifique est Fequation d’etat du fluide a temperature nulle, e’est-a-dire ici la densite volumique d’energie eo(p) de 
Fetat fundamental du systeme spatialement homogene de densite p. 

Afin d’obtenirla derivee temporelle de F operateur phase du condensat de paires de fermions, nous apportons deux 
raffinements a la theorie, une regularisation des divergences ultraviolettes et une specialisation au cas irrotationnel. 

D’une part, nous resolvons les problemes de divergence du niveau d’energie fundamental du hamiltonien de 
Landau-Khalatnikov, dus a Fenergie quantique de point zero des modes propres du systeme, en discretisant Fes- 
pace en un reseau cubique de pas Z, sous-multiple de la longueur L de la boite de quantification, beaucoup plus petit 
que la longueur d’onde typique 2nlqih des excitations du fluide peuplees thermiquement, mais beaucoup plus grand 
que la distance moyennep entre particules, 

p-'/3 « I « (36) 


ce qui est realisable a suffisamment basse temperature; ce maillage souscrit au regime de validite de Fhydrodyna- 
mique, celui d’une description spatiale a gros grains, et fournit une coupure naturelle ultraviolette dans Fespace des 
vecteurs d’onde q, en les restreignant a la premiere zone de BrillouinpiD = [-njl, njl^. Dans le hamiltonien, il faut 
alors remplacer les operateurs differentiels (gradient, divergence, laplacien) par leurs equivalents discrets, comme il 
sera implicitement fait plus bas, et il faut utiliser la densite volumique d’energie nue eo^ip), fonction du pas du re¬ 
seau /; suivant les idees de la renormalisation, Fenergie quantique de point zero des modes, formellement divergente 
lorsque I — » 0, s’ajoute alors a eo,o(p) pour reconstruire precisement la densite volumique d’energie effective ou vraie 
eoip), independante de I, et qui est celle mesuree experimentalement. 

D’autre part, nous specialisons la theorie au cas d’un operateur champ de vitesse v(r, t) irrotationnel, qu’on pent 
done ecrire comme le gradient d’un operateur champ de phase ^(r, f), canoniquement conjugue a Foperateur champ 
de densite p(r, f): 


}% S 

v(r, f) = — grad0(r, f) avec [p(r,t),^(r',t)] - 

m P 


(37) 


Ceci revient a negliger la composante transverse du champ v(r, f), comme Fa fait la reference Il22ll pour determiner 
Finteraction entre phonons et avancer dans le calcul de la viscosite. Dans le cas particulier invariant d’echelle du gaz 
de fermions unitaire, une justification en est donnee dans la reference 1351] dans le cadre de la theorie effective des 
champs. Notons en passant que la densite de quasi-particules fermioniques, exponentiellement petite en 1/T, est elle 
directement omise par la theorie hydrodynamique. 

La suite est assez classique. On part des equations du mouvement des champs en representation de Heisenberg, 
respectivement Fequation de continuite quantique et Fequation d’Euler quantique sur le potentiel (dont le gradient 
donne Fequation d’Euler quantique sur la vitesse): 


d,p + div 


1 




= 0 


= --mV - po,o(p) 


(38) 

(39) 


19. Il faut aussi que les excitations elementaires du systeme restent phononiques sur tout le domaine 2). Ceci impose hc/l < A done mcl/h > 
mc^ j/S. » 1 dans la limite de BCS, et seulement mclih > 1 ailleurs, ou c est la vitesse du son et A le gap aT = 0. La condition necessaire mclih > 1 
est plus faible que celle deja invoquee, » 1, dans le regime d'interaction forte ou dans la limite de BCS, puisque c y est d'ordre la vitesse 
de Fermi. Elle est plus contraignante dans la limite dite de CBE, ou les paires de fermions, assimilables a des bosons, forment un condensat en 
interaction faible. 
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ou {A, B} - AB + BA est Fanticommutateur de deux operateurs et 


Bo.oip) = ^eo,o(p) = e'Qoip) 


( 40 ) 


est le potentiel chimique nu de Fetat fondamental a la densite p. Ces equations derivent du hamiltonien 


^hydro ~ ^ ^ ] 


1 

-m\ ■ p\ + eofiip) 


(41) 


Les fluctuations spatiales de densite et de phase sont faibles, parce que p03/ assez grand si elles sont quantiques, 
parce que est de plus assez petit si elles sont thermiaues.l^ On pent done, comme dans |22], lineariser les equations 
du mouvement autour de la solution spatialement uniforme : 


p(r, f) = po + 6p(r, t) (42) 

0(r,f) = 0o(O + <^0(r, 0 (43) 

L’operateurpo vaut simplement NIL?, ou N est Foperateur nombre total de particules ; e’est une constante du mouve¬ 
ment. L’operateur est Foperateur phase du condensat, ici 

00 = 6o/2 (44) 


puisque Foperateur phase Oq de Fequation © prend les paires pour objets elementaires, alors que les equations 
(I38I39I) privilegient les particules fermioniques. Les fluctuations spatiales 6p et 50, d’integrale (discrete) nulle sur tout 
Fespace, se developpent seulement sur les ondes planes de vecteur d’onde q non nul, et commutent avec po- Apres 
resolution des equations linearisees sur dp et 50, on dispose de Fhabituel developpement sur les modes propres : 


-vl/2 

5p(r,0 = ^ Z 


hq 




2mco{ 


1/2 


(fiq + 


50(r, t) 




Po 


i_ Y' ( 

r3/2 Zj I 97?/; 


4/2 


qe^Z3.n» 


1 2/1^ / 


(Bq- 


Blq)e‘‘i-- 

(45) 

Blq)e“i-- 

(46) 


oil les operateurs de creation B^ et d’annihilation Bq d’un phonon de vecteur d’onde 
aux relations de commutation bosoniques, par exemple [Bq, B^,] = <5q,qs et ou Fon a 
son nue a temperature nulle 


co,o = 


|PoBo,o(Po)j‘^^ 


q et d’energie Hqcofl obeissent 
introduit Foperateur vitesse du 

(47) 


Ilreste alors a developperle secondmembre de (l39l) iusau’aFordredeuxinclus en5p et 50, a en extrairelacomposante 
de Fourier de vecteur d’onde nul, a effectuer une moyenne temporelle a gros grains pour supprimer les termes croises 
oscillants BqB_q et BlqBq, puis a reconnaitre Fidentite 


^[PBo,o(P)]^^^ = 




(48) 


20. On verifie d’abord que les fluctuations de densite sur un meme site sont faibles en valeur relative, en utilisant ED. A r = 0, <5p2)/p2 ~ 
<g: 1, compte tenu de ED et de la condition necessaire mcilh > 1 etablie dans la note 19. A T >0 s’ajoute la contribution thermique 

« ^ meme <«: 1. On verifie ensuite que les fluctuations de phase entre deux sites voisins sont faibles en valeur absolue. Pour cela, 


on constate d’abord, en utilisant ED, que P((grad(5^)^) = {mcljh)^{6p^)!. II reste ensuite a utiliser la propriete = 0(1) qui vaut 

dans toute la zone de raccordement entre les limites de CBE et de BCS, ainsi que nos estimations precedentes de pour conclure. 
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pour obtenir 


=-mo(po) - Xi \^:^j 


qe 


^I?T\D 


AN 


-Co,0 


BlK 


1 

+ — 

2 


(49) 


A cet ordre du developpement en les fluctuations spatiales, on pent regrouper dans ( l49l l la contribution de point zero 
des modes (le terme 1/2 entre parentheses) avec le potentiel chimique nu //o,o(po) pour former le potentiel chimique 
vrai jUoCpo) du fluide a temperature nulle, et identifier co,o dans le prefacteur des operateurs nombres de phonons BqBq 
avec la vitesse du son vraie a temperature nulle, co = [poA(o(po)/™]*^^- On obtient alors la limite phononique (de basse 
temperature) de la relation ([8]), et ceci sans hypothese sur la force (ou la faiblesse) des interactions. 
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